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DIFFERENTIALE DER BRECHKRAFT EINES BELIEBIGEN OPTI- 

SCHEN SYSTEMS NACH SEINEN KONSTRUKTIONSELEMENTEN, 

KRÜMMUNGEN, LINSENDICKEN, LUFTABSTANDEN UND BRE- 
CHUNGSZAHLEN 


Von J. BARTKOWSKA 
Polnische optische Werke, Warszawa 
(Eimgegangen am 27 Januar, 1961) 


Die Brechkraft eines Systems von unendlich dünnen Linsen wird nach Brechkräften von einzelnen Linsen 
differenziert. Die Brechkraft eines beliebigen optischen Systems, welches aus Plan- und Kugelflächen besteht 
wird als Funktion seiner Konstruktionselementen behandelt. Es werden einfache Formeln für erste und für 
höhere Ableitungen und Differentiale gegeben. 


1. Differentiale eines Systems von unendlich dünnen Linsen 


Die bekannte Brechkraftformel für ein beliebiges System von unendlich dünnen Linsen 
lautet 


P= X hp 


« die Brechkraft des ganzen Systems 
gy die Brechkräfte von einzelnen Linsen. 
Dasselbe Resultat wird beim Vorwärts- wie beim Rückwärtsrechnen erhalten 


P= 2 hp (1) 
دب‎ (2) 
Es wird in beiden Fällen ,1ع يم‎ s,— oe gesetzt. Im folgenden werden alle Daten aus der 
Vorwärtsrechnung bei s,— oo und =l mit dem Vorwärtspfeilchen (Abb. 1) und alle 
Daten aus der Rückwärtsrechnung (Abb. 2) mi dem Rückwärtspfeilchen bezeichnet. Wir 
wollen beweisen, dass 
UL TE (3) 
09 
Die Ableitung ist gleich dem Produkte der Höhen des Vorwärts- und des Rückwärtsstrahles. 
Alle Indices sind weggelassen, weil beide Höhen der betrachteten Linse entsprechen. Um die 
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Abb. 1. Der Vorwärtsstrahl 


Abb. 2. Der Rückwärtsstrahl 


Formel (3) zu beweisen werden wir bekannte Formel fiir die Brechkraft eines Systems 
welches aus zwei Teilen besteht, benótigen 


P = pa T Pg — epa PB (4) 
e bedeutet hier den Abstand zwischen entsprechenden Hauptpunkten beider Teile. Der 
Hauptpunktenabstand besteht aus dem Abstande d von der letzten Linse des A-Teiles bis 


zur ersten Linse des B-Teiles und beiden Hauptpunktschnittweiten (Abb. 3). Die letzten 
sind gleich den Unterschieden zwischen inneren Brennweiten und Schnittweiten 


e=d+ fa —fa ha fg —/a ha 


| Hy | Ha 
| | 
| 


Abb. 3. Zum Beweis der Brechkraftformel (5) 


h 4 bedeutet dabei die Höhe des Vorwärtsstrahles auf der letzten Fläche des 4 -Teiles, hp die 
Hohe des Rückwärtsstrahles auf der ersten (inneren) Fläche des B-Teiles. Die Produkt 
Sa h 4 und fp h B bedeuten hier die inneren Schnittweiten. Aus der Formel (4) folgt: 


p = hp Pat ha YB — APA 9p (5) 
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Das ganze optische System wird so geteilt, dass die veränderliche Linse zum B-Teile als 
erstes Glied gehört, hp bedeutet jetzt die Höhe des Strahles auf dieser Linse, hy auf der vorigen 


’ Linse (Abb. 4). In der Formel (5) sind hy, ^5 und 94 konstant und nur pg veränderlich. Der 
Differentialquotient ist gleich 


(o + oy) — p(y)  (ha—dga)ógs موقو‎ 
Oy 09 Oy 


— + — دل 
a‏ 

RS So 

` 


Abb. 4. Zum Beweis der Formel (3) 


+ 


Ds = » — dq 4 nach der bekannten Übergansformel. Die Formel (2) wird für das B-Teil ange- 
wendet und differenziert. Es folgt 


Ógg = hg 09 


woraus die Formel (3) unmittelbar folgt. 
Es wurde kein Grenzübergang für kleine 6g durchgeführt und so bleibt die Formel 


dep = hh ôy (6) 


für beliebig grosse ðw richtig. Alle höheren Ableitungen nach derselben Veränderlichen 
verschwinden. Aus der Formel (5) folgt noch 


IP _ 
Sup VWE 


Auch 
Op = — Pap dd 
bleibt richtig für beliebig grosse ðd. 


2. Die Differentiale der Brechkraft und der Brennweite nach Flächenkrümmungen 


Die Formel (3) lisst sich für ein beliebiges optisches System, welches aus Kugel- und 


Planflächen besteht, verallgemeinern als 


Ip + 
2o hhAn (7) 
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oder 


dep = hhAndo (8) 


0 = N nedeutet die Krümmung. 
r 


Auch die Formel (8) bleibt richtig für beliebig grosse 6g. Die Brechkraft einer Spiegelfläche 


ist gleich 


P 2NOsp 


nist nur dann von 1 verschieden, wenn die Spiegelung an der Hinterfläche einer Linse erfolgt. 
Die Ableitung der Brechkraft nach der Krümmung einer Spiegelfläche ist gleich 


A 
A 2T 
kad Ta 2nh h 


ó) Osp 


und das Differential 
op = — Inhhög,, 


Planflächen gibt es nur theoretisch. In der Wirklichkeit besitzt jede „Planfläche“ eine mess- 
bare Brechkraft óp oder Krümmung ðo. Für „Planflächen“ eignen sich gut die Formeln (3) 
und (6) oder (7) und (8). 


Die Ableitungen der Brechkraft und Brennweite nach Radius haben folgende Gestalt 


Ip >- An 
— = hh, 
Or r 
OF F2. 
A =—hhdn 
er F 
Die hóheren Ableitungen welche nicht mehr verschwinden, lassen sich nach gewðhnlichen 


Formeln der Differentialrechnung bestimmen. Die Näherungsformel 


DR 


r P 


ist gut für Toleranzrechnungen geeignet, 


3. Differentiale nach Linsendicken und Luftabständen 


Die veränderliche Dicke einer Linse teilt das optische System in zwei Teile. Die erste 
Fläche dieser Linse gehört noch zum A-Teile die zweite zum B-Teile. Die Brechkraftformel 
(5) hat jetzt folgende Gestalt 


- E d 
P = hg pa + ha QB — „ PAPB (9) 


Die Brechkraft des A-Teiles p a lässt sich von der Vorwärtsrechnung, des B-Teiles pp von 


X2 


der Rückwürtsrechnung, ablesen. 
Q4 na (10) 
PB = na (11) 


& und & bedeuten Winkel zwischen Strahl und Achse. Die Ableitung ist gleich 


^ 
2 — — n&& (11a) 
d,na und & liegen in demselben Raume. 
Die Formel 
dp = — naadd (11b) 


bleibt richtig für beliebig grosse ód. Für die Differentiale nach Luftabständen wird n=1 
gesetzt. 
Fiir kleine Linsendicken oder Linsenabstandinderungen haben wir niiherungsweise: 


ôF >+ 

Fe = n& aod 
Es ist leicht zu erkennen, dass 

OF 

(e7 - 

od = Fn 

© A 


4. Differentiale nach Brechungsindex 


Der Raum, welcher dem veränderlichen Brechungsindex entspricht, teilt das optische 
System i in zwei Teile A und B (Abb. 5). Wenn wir die Formel (9) nach n differenzieren (nur 


hp und hy diirfen als constant betrachtet werden), bekommen wir 


Ip - d IPA > d ونوك‎ , d: PAPB 12 
on ` g n o) xð (i ENT An? az) 
Wenn wir die Formeln (1) und (2) zum 4 und B-Teile anwenden und nach n differenzieren, 
haben wir 
Spa _ ha (13) 
on 74 
S (14) 
on TB 


Beim Rückwärtsrechnen werden Radienvorzeichen geändert. Wenn man noch die Formel 


(10) und Übergangsformel 
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Abb. 5. Die Teilung eines optischen Systems beim Differenzieren nach der Dicke oder Brechungszahl der zweiten 
1 Linse 


berücksichtigt, erhält man nach einigen Umformungen 
++ + æn 
3 hąh hgh 
c 2 B'*B - >+ 
Ha + dać (15) 
on TA TB 


Der Index 4 entspricht der ersten, und B der zweiten Fliiche der betrachteten Linse. Die 
Näherungsformel 


go [Fahr _ hohe 
3 TA TB 


+ dż | ón (16) 


gibt in den praktischen Rechnungen so gute Resultate, dass sie in genauen Abberations- 
rechnungen Anwendung findet. 


5. Höhere Differentiale 


Die zweite Ableitung nach demselben Brechungsindex durch Differenzieren der Formel 
(12) erhalten. Aus den Formeln (13) (14) und (10) (11) folgt nach einigen Umformungen 


8hp [AT ng 
r-|--à||-2-.x|— 
On? TA TB ln 
Die zweite gemischte Ableitung nach zwei benachbarten Kriimmungen findet man durch 
das Differenzieren der Formel (7). Die Flächen werden mit A und B bezeichnet. Wenn wir 


die Formel (7) für die Flåche B anwenden und hp , welches von g, abhängig ist, mit der Über- 
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gangsformel, ersetzen, bekommen wir 

Ip d 
el Agn 
Dox B | RA ae B 


Die 4-Fläche ist die letzte im A-Teile. Wenn man die Formel (12a) nach 04 differenzieren, 
haben wir 


0 Lad 
n = haan 
9 04 
Daraus folgt 
0¢ Vi | , 
00400B = — hąhgdand gn = m 


wo d und n im Raume zwischen beiden Flächen liegen. Im besonderen Falle, wenn beide 
Flächen eine in Luft befindliche Linse bilden, bekommt man die Formel 
op ++ (n—ly*d 


20268 hale n 


(17a) 
Die Formel (17) enthält keine von o4 oder og abhängige Glieder. Alle höheren Ableitungen 
nach denselben veränderlichen o4 und op verschwinden. Mittels Formeln (7) und (17a) wird 
eine genaue Formel gegeben für die Änderung der Brechkraft eines beliebigen optischen 
Systems, in welchem gleichzeitig zwei Krümmungen einer Einzellinse geändert wurden. Die 
folgende Formel bleibt richtig für beliebig grosse ðo, und 603, 


din E: 


in = 1) (hahador = hgh BOOB) JE hahp 604008 


Die Formel (17a) kann weiter differenziert werden. So findet man héhere Ableitungen nach 
Dicke und Brechungsindex derselben Linse: 


9€ ++ (n—1)? 
TL 


Ppop a 
p > + l ) 
2020897 — hahzd (1 E n? 


Wir geben noch ohne Beweis die Formel für die zweite Ableitung der Brechkraft nach Krüm- 
mung und folgenden Abstand oder Linsendicke an: 


Ip E 
309d =h gAn 
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REZIPROZITATSTHEOREM UND BABINETSCHES PRINZIP IN DER 
KIRCHHOFFSCHEN THEORIE DER BEUGUNG 


VoN A. RUBINOWICZ 
Polnische Akademie der Wissenschaften 
(Eingegangen am 23. März 1961) 


Auf Grund der Symmetrieeigenschaften der einfallenden und der gebeugten Welle gegeniiber 
der Vertauschung der Lagen des Beobachtungspunktes und der Lichtquelle sowie der Lagen 
zweier komplementärer beugender Schirme, werden das Reziprozititstheorem sowie das Babi- 
netsche Prinzip für Kirchhoffsche Wellenfunktionen nicht nur, wie dies gewöhnlich geschieht, 
für den Schatten- sondern auch für den Licht-Halbraum angegeben. Es wird jedoch gezeigt, daß 
unter Zugrundelegung der in den folgenden Überlegungen benutzten Definition der Kirchhoff- 
schen Wellenfunktionen in den beiden mathematisch ganz einwandfreien Reziprozitätstheoremen 
für den Licht- und Schatten-Halbraum auch Kirchhoffsche Wellenfunktionen auftreten können, 
die keine sinnvolle physikalische Bedeutung haben. 


Im folgenden soll gezeigt werden, wie im Rahmen der Kirchhoffschen Theorie der: 
Beugung das Reziprozitätstheorem und das Babinetsche Prinzip lauten, wen sie für alle 
Raumpunkte und nicht nur für spezielle, wie dies gewöhnlich geschiecht (vgl. z. 
B. Rubinowicz 1957), formuliert werden. Wir behandeln diese beiden Sätze gemeinsam, 
weil sie hier bei den Sprungwertproblemen der Schwingungsgleichung eine gemeinsame 
Wurzel haben. Diese wird durch die Symetrieeigenschaften der beiden Wellenbewegun- 
gen gegeben, in die eine Kirchhoffsche Wellenfunktion aufgespalten werden kann, 
nämlich durch die der direkt einfallenden und die der gebeugten Lichtwelle. Und zwar 
handelt es sich hier um zwei verschiedene Symmetrieeigenschaften der genannten. beiden 
Wellenfunktionen. Die eine tritt bei der Vertauschung der Lage von Lichtquelle und 
Beobachtungspunkt auf. Die andere gilt jedoch bei der Vertauschung von zwei kom- 
plementüren beugenden Schirmen. Diese Symmetrieeigenschaften sind so beschaffen, dab 
auch Relationen auftreten, die eigentlich ein Mittelding zwischen dem Reziprozitätstheo- 
rem und dem Babinetschen Prinzip darstellen, da sie zwei Kirchhoffsche Wellenfunk- 
tionen enthalten, die sich nicht nur durch die Vertauschung der Lage von Lichtquelle 
und Beobachtungspunkt sondern auch der Lagen von zwei komplementären beugenden 
Schirmen voneinander unterscheiden. 

Diese Sachlage ist durchaus verschieden von der, die wir bei den als Randwert- 
aufgaben im Gebiete der Schwingungsgleichung formulierten Beugungsproblemen antref- 
fen (vgl. $ 5.) Hier stellt das Reziprozitätstheorm die Symmetrieeigenschaften der gesamten 
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Wellenfunktion bei der Vertauschung der Lage der Lichtquelle und des Beobachtungs- 
punktes dar. Das Babinetsche Prinzip beschreibt jedoch die Symmetrieeigenschaften der 
gesamten Wellenfunktion bei der Vertauschung von zwei komplementiren Schirmen, 
wobei an dem einen Schirm die Randbedingung u = 0 an dem anderen aber Qu/dn = 0 
vorgeschrieben ist. 


$ 1. Symmetrieeigenschaften der einfallenden und der Beugungswelle 


Sei L eine harmonisch periodisch schwingende Lichtquelle und 5 ein allseits ins Unend- 
liche sich erstreckender, flächenhafter Schirm. In S soll sich eine, durch einen beugenden 
Rand B begrenzte, ganz im Endlichen liegende Beugungsöffnung befinden. Wir denken 
uns sie durch eine den beugenden Rand B überspannende Fläche, d. h. durch den komple- 
mentären beugenden Schirm S’ verschlossen. Durch die beiden Flächen S und 5’ wird der 
ganze unendliche Raum (im folgenden abgekürzt mit g. R. bezeichnet) in zwei Halbräume 
aufgeteilt. Der eine Halbraum, der die Lichtquelle Z enthält, soll als der Lichthalbraum 
(L. H. R.) und der andere als der Schattenhalbraum (Sch. H. R.) bezeichnet werden. Dabei 
wollen wir, der Kürze des Ausdruckes wegen, übereinkommen zu sagen, daß diese oder jene 
Beziehung im L. H. R. oder Sch. H. R. oder auch im g. R. gilt, falls dies für die dort befin- 
dlichen Beobachtungspunkte P der Fall ist. 

Nach Kirchhoff (vgl. etwa Rubinowicz 1957) wird dann im Falle eines nicht reflektieren- 
den Schirmes S die den ganzen unendlichen Raum erfüllende Wellenbewegung angenähert 
durch den Ansatz 


N | eikr I [eike eike 9 [gihr 
ur(L,P;S)= IR 7 (lu (1.1) 
5 — چ‎ — 
beschrieben. R, r und o stellen hier die Absolutbeträge der Vektoren R = LP, r = PQ, p= LQ 
dar, wobei Q die Punkte der Schirmfläche S bezeichnet, über die die Integration erstreckt 
wird. Die Normale n an den Schirm S ist hier in den Sch. H. R. gerichtet. Dabei setzen wir 
im folgenden voraus, daß die Integration ausnahmslos über den ganzen beugenden Schirm S 
erstreckt wird, also auch über seine nicht beleuchteten Partien, falls solche, etwa wenn S 
gewellt ist, auftreten. Diese Annahme stimmt nicht ganz mit der ursprünglichen Kirch- 
hoffschen Voraussetzung, daß die Integration nur über die beleuchteten Teile des Schirmes 
5 auszuführen ist. Sie trägt aber wesentlich zur Einfachheit der Formulierung der weiter unten 
anzugebenden Sätze bei, ist aber gleichzeitig eine Quelle von Paradoxien, die wir in $5 
besprechen. 


Die durch den Kirchhoffschen Ansatz (1.1) vom Fresnelschen Standpunkte beschriebene 
Wellenbewegung läßt sich vom Youngschen Gesichtspunkte aus in der Gestalt 


ug(L, P; S) = ug(L, P; S) + ug(L, P; S) (1.2) 

darstellen (Rubinowiez 1917, 1957). , 
ug(L, P; S) ist hier die einfallende Lichtwelle. Sie ist im Lichtraum, d. h. in allen Beo- 
bachtungspunkten P, von denen aus man die Lichtquele L sehen kann, gleich e'**/R und 


verschwindet im Schattenraum, d. h. in allen übrigen Raumpunkten P. Sie beschreibt also 
die Lichtverteilung gemäß den Gesetzen der geometrischen Optik. 
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ug(L, P; S) stellt die Beugungswelle dar, die durch das über den beugenden Rand B 
erstreckte Integral 


l | eee r(p x ds) 
de, PoS) = 
up( 5) a A eS (1,3) 
B 


gegeben wird. ds bedeutet hier ein Bogenelement des beugenden Randes B. Der Richtungs- 
sinn von ds wird dabei durch die Rechtsschraubenregel mit der Normalen n’ an die die 
Beugungsöffnung verschließende Fläche S’ gekoppelt. Dabei wird angenommen, daß die 
Normale n’ in den L. H. R. der Lichtquelle Z hinein gerichtet ist. 

Wir betrachten zunächst die Symmetrieeigenschaften der direkt einfallenden und der 
Beugungswelle bei der Vertauschung der Lage der Lichtquelle Z mit der des Beobachtungs- 
punktes P, die wir vor allem für die Begründung des Reziprozitätstheorems benötigen. 

Zunächst ist es klar, daß die einfallende Lichtwelle ug(L, P; S) das Reziprozitätstheorem 


ug(L, P; S) = ug(P, L; 5) im g. R. (1.4) 


erfüllt. Befindet sich nämlich die Lichtquelle Z in Sicht des Beobachtungspunktes P, so liegt 

auch P in Sicht von L und sowohl die linke als auch die rechte Seite der Beziehung (1.4) 

werden durch e®®/R gegeben. Ist aber L von P aus nicht sichtbar, so kann auch P von L 

aus nicht gesehen werden. Die beiden in (1.4) auftretenden Funktionen verschwinden sodann. 
Für die Beugunswelle ug(L, P; S) (1.3) gilt das Reziprozitätstheorem 


Bla OS) nU LS) im sch. HH. R. (1.5) 


Das Vorzeichen der Beugungswelle (1.3) bleibt nämlich bei der Vertauschung von L und P 
ungeändert, weil ja dabei nicht nur die Rollen von r und م‎ vertauscht werden, sondern auch 
der Richtungssinn von ds umgekehrt werden muß. Wir müssen ja doch den Richtungssinn 
der Normalen n’ in den entgegengesetzten abändern, weil er ja doch stets in den L. H. R. 
der Lichtquelle Z hinein gerichtet sein muß und der Richtungssinn von ds mit dem von n' 
durch die Rechtsschraubenregel verkniipft ist. 

Liegt jedoch P im L. H. R., so behalten n’ und ds bei der Vertauschung von L und P 
ihren Richtungssinn bei und die Vertauschung von r und p bewirkt eine Vorzeicheninderung 
der Beugungswelle, so daß 


UAL, P; 3) == uz, 255) im L. H.R. (1.6) 


Die Beugungswelle (1.3) ist somit im Sch. H. R. eine symmetrische und im L. H. R. 
eine antisymmetrische Funktion bei der Vertauschung der Lage des Beobachtungspunkies P 
und der Lichtquelle L. 

Wir müssen nun, vor allem für die Ableitung des Babinetschen Prinzips, angeben, wie 
sich ug(L, P; S) und ug(L, P; S) ändern, wenn wir von dem Kirchhoffschen Beugungsproblem 
für den Schirm Szu dem komplementären Beugungsproblem für den Schirm 5" übergehen. 

Für die einfallende Lichtwelle gilt dann 


ug(L, P5 S) UNL P; 5") —e*RIR im L. H. R. (1.7) 
und 


lL, ONS mL. ES e^RIR. im Sch. H. R. (1.8) 
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Im L. H. R. wird ja die einfallende Lichtwelle, unabhiingig von der Gestalt des beugen- 
den Schirmes S, durch die unbehinderte Wellenbewegung des von der Lichtquelle L ausge- 
strahlten Lichtes gegeben. Im Sch. H. R. ist jedoch die Summe der einfallenden Lichi- 
wellen für die beiden komplementären Schirme in allen Punkten dieses Halbraumes gleich 
der unbehinderten Wellenbewegung der Lichtquelle. 


Für die Beugungswelle (1.3) gilt jedoch die Beziehung 
u,(L, P; S) = — ug(L, P;S) im g. R. (1.9) 


Beim Übergang von der Beugungswelle ug(L, P; S) für den Schirm 5 zu der Beugungswelle 
ug(L, P; S’) für den komplementären Schirm S müssen wir nämlich, unabhängig von der 
Lage des Beobachtungspunktes P den Umlaufssinn auf dem beugenden Rande B um- 
kehren. Während ja der Richtungssinn von ds im Falle von ug(L, P; S) mit der nach dem 
L. H. R. der Lichtquelle Z hin gerichteten Normalen n’ auf 5’ durch die Rechtsschrauben- 
regel gekoppelt war, ist er im Falle von ug(L, P; S’) mit der nach dem Sch. H. R. der Licht- 
quelle L hin weisenden, auf dem Schirm S befindlichen Normalen durch die gleiche Regel 
verknüpft. 


§ 2. Reziprozitätstheorem 


Zunichst beschiftigen wir uns mit dem Reziprozititstheorem. In dem von v. Laue 
(1928) betrachteten Falle, wo der Beobachtungspunkt P sich im Sch. H. R. befindet, folgt 
aus (1.2) sowie aus der Tatsache, daß gemäß (1.4) und (1.5) hier die direkt einfallende und 
die Beugungswelle das Reziprozitätstheorem erfüllen, daß dieses Theorem auch für die 
Kirchhoffsche Lósung (1.1) hier gilt: 


ur(L, P; S) = ug(P, L; S) im Sch. H. R. (2.1) 
Im Lichthalbraum erhalten wir jedoch durch Addition von (1.4) und (1.6) 
ug(L, P; S) + ug(L, P; S) = ug(P, L; S) — ug(P, LS S) im L. H. R. (2.2) 


Mit Rücksicht auf (1.7) kann jedoch auf der rechten Seite von (2.2) die Wellenbewegung 
ug(P, L; S) durch ug(P, L; S") ersetzt werden. Da gemäß (1.9) ug(L, P; S) beim Ersatz | 
von S durch S' jm g. R. sein Vorzeichen ändert, so kann in (2.2) ug(P, L; S) durch 
—u,(P, L; S') ersetzt werden. Die rechte Seite von (2.2) ist somit gleich ug(P, L; S") 
+u,(P, L; 5’), so daß mit Rücksicht auf (1.2) die Beziehung (2.2) in 


Ug(L, P; S) = ug(P, L; S" im L. H. R. (2:53) 
übergeht. 

Liegt also der Beobachtungspunkt P im Sch. H. R., so ändert sich gemäß (2.1) nicht 
die Wellenbewegung bei der Vertauschung der Lage des Licht- und Beobachtungspunktes, 
L bzw. P. Befindet sich jedoch der Beobachtungspunkt P im L. H. R., so erhalten wir mit 
Rücksicht auf (2.3) bei der Vertauschung der Lage von L und P die Wellenbewegung für 
den komplementären beugenden Schirm S’. 

Da in der Beziehung (2.3) nicht nur die beiden Punkte L und P, sondern auch die beiden 
Schirme S und S’ miteinander vertauscht werden, so stellt diese Beziehung eigentlich ein 
Mittelding zwischen dem Reziprozitütstheorem und dem Babinetschen Prinzip dar. 
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$ 3. Babinetsches Prinzip 


Das Babinetsche Prinzip der Kirchhoffschen Beugungstheorie wird gewöhnlich für den 
Fall ausgesprochen (vgl. etwa Rubinowicz 1957), wo der Beobachtungspunkt P sich im 


Sch. H. R. befindet. Es lautet dann 
Ug(L, P; S) + ug(L, P; S") = e*F/R im Sch. H. R. (25 


Um den Beweis von (3.1) zu erhalten, braucht man bloß zum Babinetschen Prinzip (1.8) 
für die einfallende Lichtwelle die Summe der Beugungswellen ug(L, P; AZYL LORNA 
zu addieren, die nach (1.9) verschwindet. 


Im L. H. R. erhält man hingegen für die in (3.1) links auftretende Summe 
UL, P; 8) + uz(L, P; P im L. H. R. (3.2) 


Die Summe der links in (3.2) auftretenden beiden einfallenden Wellen ergibt ja gemäß 
(1.7) die rechte Seite dieser Beziehung, während die Summe der entsprechenden beiden 
Beugungswellen mit Rücksicht auf (1.9) verschwindet. 

Während also die Summe der beiden Wellenbewegungen, die durch den Kirchhoffschen 
Ansatz (1.1) für zwei komplementäre Schirme gegeben werden, gemäß (3.1) im Sch. H. R. 
die einfache einfallende Welle ergibt, erhält man für diese Summe im L. H. R. laut (3.2) 
den doppelten Wert. 

Im Falle der Fraunhoferschen Beugungserscheinungen sind die rechten Seiten von 
(3.1) und (3.2) gleich Null zu setzen, so daß die Intensitäten der Beugungserscheinungen 
sowohl im Sch. H. R. als auch im L. H. R. jeweils für zwei komplementäre Schirme S und S$”: 
einander gleich sind. Dies ist vom Youngschen Standpunkt selbstverständlich, weil ja im 
Fraunhoferschen Grenzfalle die Beugungserscheinungen in jedem der beiden Halbräume 
allein durch die Beugungswelle ug(L, P; S) (1.3) verursacht werden. Die einfallende Licht- 
welle ug(L, P; S) geht nämlich in diesem Grenzfalle in eine ebene Welle von einer ganz 
exakt festgelegten Fortschreitungsrichtung über, die in den Sch. H. R. hirein gerichtet ist. 
Sie kann sich also bei der Beobachtung des Fernfeldes der Fraunhoferschen Beugungser- 
scheinungen nur in einem bestimmten Punkte der Beobachtungsebene im Sch. H. R. be- 
merkbar machen, wo ihr die Aufgabe zufällt eine von der Schattengrenze der Fresnelschen 
Beugungswelle stammende Unstetigkeit zu glitten (vgl. Rubinowicz 1954, 1957). Im Licht- 
halbraume kommt sie überhaupt nicht zur Geltung, da sie ja hier in der entgegengesetzten 
Richtung fortschreitet wie die vom beugenden Rande ausgehenden, zur Beobachtung ge- 
langenden Beugungswellen. 

Allerdings muß bemerkt werden, daD nach der Kirchhoffschen Theorie die Fraun- 
hoferschen Beugungserscheinungen im L. H. R. viel lichtschwicher sein sollten als im 
Sch. H. R. weil sie ja, wie oben bemerkt wurde, in den beiden Halbriumen eigentlich nur 
durch die Beugungswelle verursacht werden und diese in der Nihe der Schattengrenze also 
im Sch. H. R. die größten Amplituden hat. 

Mit Hilfe der oben angegebenen Symmetrieeigenschaften der einfallenden und der 
Beugungswelle lassen sich auch noch andere Beziehungen von der Art der oben angegebenen 
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ableiten. Als Beispiele führen wir an 
ug(L, P; S) + ug(P, L; S) = 2e"R/R im L. H. R., (3.3) 
ug(L, P; S) — ug(P, L; S) = 2ug(L, P; S) = — 2ug(P, L; S) im LHR. (34) 
ug(L, P5 S) — ug(L, P; S) = 2ug(L, P; S) = — 2ug(P, L; $) im L.H.R. (3.5) 


Dafür, daß man die Beziehung (3.3) eventuell als eine Abart des Reziprozitätstheorems 
für den L. H. R. ansieht, spricht die Tatsache, daß in ihr zwei Kirchhoffsche Wellenfunktio- 
nen auftreten, in denen die Rollen der beiden Punkte P und Z vertauscht sind. Dagegen 
spricht jedoch entschieden der Umstand, daß in ihr die einfallende Lichtwelle vorhanden 
ist, deren Auftreten in einem Reziprozitätstheorem man wohl nicht rechtfertigen kann. 
Die Beziehungen (3.4) und (3.5) bieten insofern kein physikalisches Interesse, als in ihnen 
die Beugungswelle ug auftritt, die durch den reichlich komplizierten Ausdruck (1.3) gegeben 
wird. 


§ 4. Andere Ableitung der obigen Beziehungen 


Es sei noch bemerkt, daß alle oben angegebenen Beziehungen auch unmittelbar aus der 
Darstellung (1.1) der Kirchhoffschen Wellenbewegung erhalten werden können. Um sie ab- 
zuleiten muß man nur beachten, daß bei der Vertauschung von L und P 

1) sich e®®/R nicht ändert, da ja R=LP ist, 

2) das in (1.1) auftretende Integral sein Vorzeichen ändert, wenn P im L. H. R. sich 
befindet und es beibehält, wenn P im Sch. H.R. gelegen ist. 

Bei der Vertauschung von L und P müssen ja im Integral in (1.1) die Entfernungen r und 
م‎ miteinander vertauscht werden, was eine Vorzeichenänderung dieses Integrales zur Folge 
hat. Befindet sich P im L. H. R., so bleibt dabei die Normalenrichtung auf S erhalten. Im 
Sch. H. R. muß sie jedoch durch die entgegengesetzte ersetzt werden, da sie ja in (1.1) stets 
nach dem Sch. H. R. der Lichtquelle hin gerichtet sein muß. Die Umkehrung des Richtungs- 
sinnes der Normalen bewirkt aber eine weitere Vorzeichenänderung des in (1.1) auftretenden 
Integrals, so daß im Ganzen dieses Integral sein Vorzeichen beibehält. 

Überdies muß man noch berücksichtigen, daß auf Grund des Helmholtz-Huygensschen 
Prinzips für den Außenraum (vgl. etwa Rubinowiez 1957, 1,3.HHa) das Integral in (1.1), 
wenn wir es über den Schirm S und zugleich auch über den komplementären Schirm S’ 
erstrecken, im Sch. H. R. (mit Rücksicht auf die Richtung der Normalen auf S und S') die 
negative einfallende Welle eRIR ergibt, im L. H. R. jedoch verschwindet. 


$ 5. Vergleich mit den analogen Beziehungen für die Randwertprobleme 


Bezeichnen wir mit ug, (L, P; S) bzw. ug. (L, P; S) die Lösung einer Randwertaufgabe 
im Gebiete der Schwingungsgleichung, die auf dem Schirme S die Randbedingung dup, [8n 
=0 bzw. ug —0 erfüllt. Das Reziprozitätstheorem für solche Randwertaufgaben bringt 
dann die Tatsache zum Ausdruck, daf die gesamte Wellenfunktion symmetrisch bei der 
Vertauschung der Lage der Lichtquelle L und des Beobachtungspunktes P ist. Es lautet 


daher bekanntlich 


Up, (L, P; S) = up, (P, L; S) im g. R. (5.1) 


Dabei kann die Gestalt des Schirmes S und die Lage des Beobachtungspunktes P eine ganz 
beliebige sein. Im Falle flåchenhafter Schirme hat somit das Reziprozititstheorem für Rand- 
wertaufgaben die gleiche Gestalt unabhingig davon, ob der Beobachtungspunkt P sich im 
Sch. H. R. oder im L. H. R. befindet. Die Gestalt des Reziprozitätstheorems (5.1) stimmt 
dabei mit dem gleichen Theorem (2.1) für die Kirchhoffsche Lösung ug (1.1) im Sch. H. R. 
überein. Im L. H. R. gilt jedoch für ug die Beziehung (2.3), die einen Zusammenhang zwi- 
schen den Kirchhoffschen Lösungen (1.1) für die beiden komplementären Schirme S und S" 
darstellt. 

Für Randwertaufgaben kann man jedoch das Babinetsche Prinzip nur im Falle ebener 
Schirme formulieren. Es stellt die Symmetrieeigenschaft der gesamten Wellenfunktion bei 
der Vertauschung von zwei komplementiren Schirmen S und S’ dar, wobei angenommen 


wird, daß an dem einen die Randbedingung dup, [9n = 0 an dem anderen jedoch up_ = 0 
erfüllt ist. Es lautet (vgl. Sommerfeld 1950, 1959): 

VANL PS) u (LP: Sys eR im Sch. H. R. (5:2) 
und 

BRI D sp 1 DCE Rm L HR, (5.3) 


Dabei bedeutet K' die Entfernung des Beobachtungspunktes P von der an der Schirmebene 
gespiegelten Lichtquelle. In allen Beziehungen (5.1), (5.2) und (5.3) gelten stets gleichzeitig, 
entweder alle oberen oder alle unteren Vorzeichen. 

Das Babinetsche Prinzip (5.2) für die Randwertaufgabe im Sch. H. R. ist sehr ühnlich 
dem gleichen Prinzip (3.1) für die Kirchhoffsche Wellenfunktion in dem gleichen Raum- 
gebiete, nur mit dem Unterschied, daß die beiden in (5.2) auftretenden, zu zwei komple- 
mentären Schirmen gehörigen Wellenfunktionen auf diesen Schirmen verschiedene Rand- 
bedingungen erfüllen. 

Im L. H. R. besteht hingegen ein betrüchtlicher Unterschied zwischen den Babinetschen 
Prinzipen (5.3) und (3.2). Im Babinetschen Prinzip (5.3) für die Randwertaufgabe tritt die 
Differenz der beiden Wellenfunktionen auf, die die Beugung an den beiden komplementären, 
jedoch verschiedenen Randbedingungen entsprechenden Schirmen beschreiben. Da aber 
diese beiden Wellenfunktionen die gleiche einfallende Lichtwelle haben, entfällt diese daher 
in (5.3) und es tritt hier nur die reflektierte Lichtwelle e*?/R' auf. Im Falle der Kirchhoff- 
schen Wellenfunktion ist kein reflektiertes Licht vorhanden und es kann sich daher im Ba- 
binetschen Prinzip nicht bemerkbar machen. Da aber (3.2) linkerhand die Summe der den 
beiden komplementären Schirmen S und S’ entsprechenden Wellenbewegungen enthält, 
muß hier mit Rücksicht darauf, daß wir uns im L. H. R. befinden, nun rechterhand die 
doppelte einfallende Lichtwelle auftreten. Die Differenz der den beiden komplementären 
Schirmen entsprechenden Kirchhoffschen Wellenbewegungen ergibt jedoch im L. H. R., 
wie die in § 3 zuletzt angegebene Beziehung (3.5) zeigt, die doppelte Beugungswelle und 
bietet daher physikalisch kein besonderes Interesse. 
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§ 6. Physikalische Grenzen für die Anwendbarkeit der bewiesenen Reziprozitätstheoreme 


Denken wir uns einen flächenhaften Schirm S gegeben, so kann in der Kirchhoffschen 
Beugungstheorie ein zu ihm komplementärer Schirm S’ einen beliebigen Verlauf haben, 
falls er nur mit dem Schirme S die gleiche Randkurve als beugenden Rand B besitzt. Bei 
gegebener Lage und Gestalt des Schirmes S kann somit jeder Beobachtungspunkt P, je nach 
dem Verlauf des komplementären Schirmes S’ in dem Sch. H. R. oder im L. H. R. der 
ursprünglich gegebenen Lichtquelle Z liegen. In den beiden angegebenen Fällen erfüllt 
jedoch die Kirchhoffsche Wellenfunktion ug (1.1) ein verschiedenes Reziprozitätstheorem 
(2.1) bzw. (2.3). 

Es ist vom mathematischen Standpunkte nicht zu bezweifeln, daß stets das Reziprozitäts- 
theorem (2.1) bzw. (2.3) anwendbar ist, sobald der Beobachtungspunkt P im Sch. H. R. 
bzw. im L. H. R. liegt, falls nur der Kirchhoffschen Wellenfunktion uz die bei allen unseren 
Überlegungen vorausgesetzte, durch (1.1) gegebene Definition zu Grunde gelegt wird, also 
in dem Integral in (1.1) die Integration ohne Rücksicht auf die Gestalt des Schirmes aus- 
nahmslos über den ganzen Schirm S oder S’ durchgeführt wird und die Normale n in den 
Sch. H. R. der betreffenden Lichtquelle hinein gerichtet ist. 

Bei der Anwendung des Reziprozitätstheorems (2.1) oder (2.3) wird man selbstver- 
ständlich stets voraussetzen, daß die ursprünglich gegebene Kirchhoffsche Wellenfunktion 
ug(L, P; S) einen physikalisch realisierbaren Wellenvorgang darstellt. Es entsteht aber die 
Frage, ob die in (2.1) bzw. (2.3) auftretende, durch Vertauschung von L und P bzw. auch 
von S und S” entstehende Funktion ug(P, L; S) bzw. ug(P, L; 5’) ebenfalls stets physika- 
lisch sinnvoll interpretiert werden kann. Um zunächst im Falle (2.1) an einem Beispiel zu 
zeigen, das dies nicht immer zutreffen muß, nehmen wir an, daß der Schirm S in einer 
Ebene E liegt, die den ganzen unendlichen Raum in zwei Halbräume H, bzw. Hç, teilt, 
von denen H, die ursprünglich gegebene Lichtquelle Z enthält. Dabei möge fürs Erste 
vorausgesetzt werden, daß der ganze komplementäre Schirm S’ in dem Halbraum H, ver- 
läuft. Der Beobachtungspunkt P möge nun in dem von 5' und der Ebene E eingeschlossenen 
Raume sich befinden, den wir etwa mit R, bezeichnen wollen (Abb. la). Da der Beobachtungs- 
punkt P im Sch. H. R. der ursprünglich gegebenen Lichtquelle Z liegt, so ist hier das Rezi- 
prozitätstheorem (2.1) anzuwenden. Tut man dies (Abb. 1b), so gelangt die Lichtquelle in 
den Raum Ro, also in den Sch. H. R. der ursprünglichen Lichtquelle. Für die Kirchhoffsche 
Wellenfunktion ug(P, L; S) einer so gelegenen Lichtquelle ist nach (1.1) eine in den L. H. R. 
der ursprünglichen Lichtquelle weisende, in der Abb. 1b mit n bezeichnete Normalenrich- 
tung auf dem Schirme S zu verwenden. Eine solche Kirchhoffsche Wellenfunktion ur(P, L; 5) 
entspricht aber keinem physikalisch sinnvollen Beugungsproblem. In der neuen Lage be- 
leuchtet nämlich die Lichtquelle den Schirm S von der gleichen Seite, wie die ursprünglich 
gegebene Lichtquelle. Soll also ug(P, L; S) einen physikalisch interpretierbaren Beugungs- 
vorgang am Schirm S darstellen, so muß man im Integral über den Schirm S die in der 
Abb. Ib mit n’ bezeichnete Normalenrichtung voraussetzen. Sie ist in den Sch. H. R. der 
ursprünglichen Lichtquelle hinein gerichtet, besitzt also einen zur Normalen n entgegen- 
gesetzten Richtungssinn, Auch als Beugungsvorgang an dem komplementären Schirme S’ 
läßt sich die in Rede stehende Wellenfunktion nicht deuten. Bei der vorausgesetzten Lage 
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Abb. 5a. ug(P, L; S) Abb. 5b. ur (P, 17357) 


Physikalische Deutung der in das Reziprozitütstheorem eingehenden Kirchhoffschen Wellenfunktionen. 
In den obigen Abbildungen werden die Schirme, über die mit Rücksicht auf die Definition (1.1) der betreffenden 
Funktion ug zu integrieren ist, durch stark ausgezogene Linien gekennzeichnet. n bezeichnet die Normalen- 
richtung, die entsprechend der Definition der betreffenden Funktion ug zu verwenden ist, n’ die Normale bei 
der die betreffende Funktion ug einen physikalisch realisierbaren Beugungsvorgang darstellen würde. Abb. 1 
und Abb. 2 beziehen sich auf das Reziprozitätstheorem (2.1), die übrigen Abbildungen auf das Reziprozitäts- 
theorem (2.3). In dem Falle der Abb. 1 bzw. Abb. 5 ist das Reziprozitätstheorem (2.1) bzw. (2.3) nicht anwendbar. 
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des Beobachtungspunktes P und des Schirmes S’ ist somit eine Anwendung des Rezipro 


zitätstheorems (2.1) physikalisch nicht zu rechtfertigen. 


Bei der angegebenen Lage des Beobachtungspunktes P, kann man jedoch, wie die 
Abb. 2a und 2b zeigen, das Reziprozitätstheorem (2.3) anwenden, wenn man nämlich vor- 
.. 9 all 5 . 9 7 „13 
aussetzt, daß der komplementäre Schirm 5’ ganz in dem Halbraume Hs, verläuft. 


Liegt der Beobachtungspunkt P im Halbraume H;s,, , so kann man, wie aus den Abb. 3a 
und 3b bzw. 4a und 4b zu entnehmen ist, sowohl das Reziprozitätstheorem (2.1) als auch 
(2.3) benützen, sobald sich nämlich der beugende Schirm S' im Halbraume H, bzw. Hsen 


befindet. 

Man kann auch Beispiele für den Fall angeben, wo bei der Anwendung des Reziprozitäts- 
theorems (2.3) bei der Vertauschung von L und P sowie S und 5’ eine Wellenfunktion 
ug (P, L; S’) entsteht, die physikalisch nicht als ein Beugungsvorgang interpretiert werden 
kann. Um ein solches Beispiel zu erhalten, setzen wir voraus, daß der komplementäre Schirm 
S' in einer Ebene E liegt und der ursprüngliche Schirm S in dem Halbraume Hs., enthalten 
ist. Der Beobachtungspunkt P möge sich in dem von der Ebene E und dem Schirm S ein- 
geschlossenen Raume befinden. Wir haben es hier somit sozusagen mit dem Falle der Abb. 
4a und 4b zu tun, in dem die Rollen der beiden Schirme S und 5’ vertauscht wurden. Wie aus 
den Abb. 4a und 4b zu entnehmen ist, hat man in der durch die entsprechenden Vertau- 
schungen entstehenden Funktion ux(P, L; 5’) in dem in (1.1) enthaltenen Integrale die Nor- 
malenrichtung n zu verwenden, während man die Normale n' benutzen müßte um einen phy- 
sikalisch deutbaren Wellenvorgang zu erhalten. 


Für die oben durchgeführten Überlegungen ist es selbstverständlich nicht wesentlich 
anzunehmen, daß der Schirm S oder 5’ eben ist. Es handelt sich hier offenbar um gewisse 
analysis situs Eigenschaften dieser Schirme, die bei der Verwendung ebener Schirme sehr 
anschaulich zu Tage treten und leicht beschrieben werden können. Die Entstehung von 
physikalisch nicht deutbaren Wellenfunktionen bei der Anwendung des Reziprozitäts- 
theorems (2.1) oder (2.3) ist somit eine auch bei nicht ebenen Schirmen S oder S’ zu erwar- 
tende Erscheinung. 


Zum Abschluß unserer Überlegungen sei noch bemerkt, daß bei nicht ebenen Schirmen 
je nach der Lage der Lichtquelle Z auch Fälle auftreten können, wo nur gewisse Teile des 
beugenden Schirmes beleuchtet sind oder auch beleuchtete Partien zu beiden Seiten eines 
beugenden Schirmes vorhanden sind. In einem solchen Falle kann man die Integration nicht 
über den ganzen Schirm erstrecken, wie wir dies bei der Definition der Kirchhoffschen Wellen- 
funktion (1.1) gefordert haben. Trotzdem könnte man wohl auch auf solche Fälle sowohl das 
Reziprozitätstheorem als auch das Babinetsche Prinzip verallgemeinern. Wir unterlassen es 
jedoch das hier zu tun, weil diese Fälle nicht sehr interessant sind und bereits an der Grenze 
der Anwendbarkeit der Kirchhoffschen Beugungstheorie liegen. In dieser Theorie gelangt 
man nämlich zu der paradoxen Schlußfolgerung, daß die in ihr auftretenden Schirme für 
die Beugungswelle vollkommen durchsichtig sind. Da ein solches Hindurchtreten der 
Beugungswelle durch einen Kirchhoffschen Schirm stets stattfinden muß, falls er nicht eben 
ist, so läßt sich diese Theorie streng genommen nur für ebene Schirme einwandfrei durch- 
führen (vgl. Rubinowiez 1957). In keinem Widerspruch steht dies zur Tatsache, daß bei 
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vorsichtiger Anwendung der Kirchhoffschen Beugungstheorie man die Beugungserschei- 
nungen auch bei nicht ebenen Schirmen behandeln kann, also auch in solchen Fallen, wo 
der beugende Rand durch eine riumliche Kurve gegeben wird. 
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ON THE MICROSCOPIC THEORY OF GIANT DIPOLE RESONANCE 


BY IW CZYŻ 
Institute of Nuclear Physics, Krakéw, Bronowice 
(Received April 4, 1961) 


The response formulation of the dipole y-ray absorbtion by closed shell nuclei is given. 
A qualitative description of the giant resonance structure including its width is proposed. 


1. Introduction 


It has been known for a long time that the so-called Giant Dipole Resonance observed 
in experiments of X-ray absorbtion by nuclei has an electric dipole nature. Several expla- 
nations accepting the collective semiclassical character of the absorbtion were proposed 
(Migdal 1944, Goldhaber, Teller 1948, Steinwedel, Jensen 1950). Later, however, the funda- 
mental role of the shell structure in nuclei became more and more evident and consequently 
a theory of the electric dipole absorbtion was given in which the giant resonance was built 
from the El single particle transitions. (Wilkinson 1956, Rand 1957). This theory was already 
a sort of “microscopic” description of the giant resonance in contradistinction to the “macros- 
copic” collective descriptions (Goldhaber, Teller 1948, Steinwedel, Jensen 1950) which did 
not take into account any internal structure of nuclei. It was soon realized that these two 
pictures were essentially equivalent (Brink 1957). There was, however, an important draw- 
back to the “microscopic” theory: it gave much too low energy of the giant resonance max- 
imum. This unsatisfactory situation recently gave an impulse to improve the shell-mo- 
del “microscopic” approach by including interactions between shell-model particles and 
holes and indeed considerable successes were achieved (Elliot, Flowers 1957, Brown, 
Bolstreli 1959, Mottelson 1960, Brown, Castillejo, Evans 1961). In these papers the explicit 
calculations of the resonance energy were made for 150 and *?Ca nuclei. 

In the present paper we try to develop a general “microscopic” approach to the pro- 
blem, based on the field theoretic formalism which has recently been applied by many authors 
to the diverse many-particle systems (see e. g. Landau 1957, Galitski, Migdal 1958, Ga- 
litski 1958, Gorkov 1958, DuBois 1959, Belyaev 1960, Gottfried, Pióman 1960 and many 


others). 
(737) 


2. Formulation of the problem 


We shall formulate our calculations in terms of the response of the nucleus to a time 
dependent external electromagnetic field". In first order (with respect to the electromagnetic 
interaction) the response is defined as follows: 


Rw) = X (w — E) 9,|H' ©, (2.1) 


where H' is the interaction between the external electromagnetic field and the nucleus, w — 
energy of the electromagnetic field and Ø,, ©, ground and excited states of the nucleus 
respectively. 

Let us define 


S(t) = <Po|T(H * (0) (0)) |,» (2.2) 


where 7(...) is Wick’s time ordering operator and the time dependence is that of the Heisen- 
berg representation. In a standard way (Goldstone, Gottfried 1959) we get 


Tm KS, H] _ ع‎ n 23 
pi) Xi w — E,4- ie w+ E, — ie CA 
W 
where S(@) = — i | dte" S(t). Hence 
żę, 1 s 
R(v) = — E Im S(w) (2.4) 


The S(w) function can provide us with all the information we need about the giant resonance: 
its poles give the resonance energies, residues at these poles being proportional to the corres- 
ponding transition probabilities. Besides, the widths of the absorbtion lines can, in principle, 
be evaluated from the analytic continuation of S(w) into one of its Riemann sheets in the 
complex w-plane (see e. g. Gottfried, Pičman 1960). 

Our calculations do not start from the “first principles” i. e. from a system of free nucle- | 
ons interacting through the real nucleon-nucleon potential. We rather take as our nonin- 
teracting system a finite nucleus with a shell structure and we let the shell model nucleons 
(which are already a kind of quasiparticles) interact through a certain effective potential. 
Consequently we expand our field operators in the shell model wave functions and satisfy 
in this way the proper boundary conditions for finite nuclei calculations. 


In the El approximation the electromagnetic interaction is: 


eh? 


/ 1 1 1 l 4 
H = Ime. dx [w30 4- Ts) 7 WY) — z Uv). (+ dy | -e (2.5) 


. . L = . . L 0 
where 73 is the standard isospin operator, é the polarization vector of the incident photon 
beam. i dx stands for integration over the space and summation over the spin and isospin 


1 The response formulation has been used by many authors in a variety 


of problems, see for example 
Gottfried, Pičman (1960) for further references. 
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variables. y is the nucleon field operator. Only the 73-containing part of (2.5) can produce 
dipole vibrations. 
We shall use the expansion 


Nast > A gP a» y* = 21 aż Pa (2.6) 
a a 


where g, belong to a set of the shell model wave function specified by the nature of our pro- 
blem. For example in the simplest situation where the giant resonance is built from the single 
particle transitions between the two neighbouring shells we limit ourselves to g, belonging 
to these two shells, stands for all the angular momentum, spin and isospin labels. We 
introduce (2.6) into the t3-part of (2.5), e being taken along the z-axis, and get: 


H' = — ieħ >) (Ep — E Jat azz(x8) (2:0) 
«b 


where E, is the energy of the Pa state and 
* 
z(o40) = | dixo rss, 
Hence the S(ż) function is: 


S(t) = eh? A B,|T(az(t)az ) (0)a;(0)) | Ø9 > x 


x (Eg — EIKE, — E,)z(x8)z(y0) (2.8) 


So, because of (2.4) and (2.8), we have reduced our problem to the investigation of the 
Fourier transform of the two-particle Green function: 


+00 
K(æpyð; o) = —i | dte**K(afyó; t) 
K(abyć; t) = <G,|T(a,(t)aż (t)a (0)a4(0)) |,» (2.9) 


We shall use the Feynman diagrams in analyzing (2.9). Our noninteracting system is 
the ladder of the shell model states filled up to a certain level «z. Our interaction Hamil- 


tonian is 


= + >] KapB|v|ðy Dag aj aa, (2.10) 


aßyö 


where 
<aplvlöy> = il d*xyd*x,9, (x3) gg) va F Xo) P(X) Pa (Ko) 
Besides, we need the Fourier transform of the bare single nucleon oo: 


0. (æ) 0 (a) 


: (2.11) 
—Btiy ل‎ 2] 


G,(4, E) = 


where 


0, ; 
9. (æ) = "n Be = E 
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We can get some idea how this formalism works calculating R(w) for a noninteracting 
system. The only graph (the *bubble") which contributes to K(w) in this approximation 


W 


I ofe E ,« 


Fig. 1 


is shown in Fig. 1. The Fourier transform of the K function is 


Ky(æpyð: w) = 0,, 0, Ky(xp: ©) 


and 
+00 
1 S 
Ky(«P ; w) = Im | aes. E+ 0)Gy(a, E) 

E 09,(0.(x) .— 0 (B)0. (æ) (2.12) 

` w—Ep+Ex+tin w+ Ek — Eg— in m 
From (2.4), (2.8) and (2.12) we get: 

Ry(w) = eh? >; 6,(B8)8_(a) (E; — E,)*z*(aB)ó(om — E,+ E,) (2.13) 


ap 


which obviously reproduces Wilkinson’s description of the giant resonance (compare Wil- 
kinson 1956). 
It would seem natural to build the dipole vibrations from the single particle excitations 


shown in Fig. 1. Therefore the function Ky(a8: w) (the particle-hole propagator) will 
play a fundamental role in our approach. 


3. A qualitative description of the giant dipole resonance for spherical nuclei 


(a) The fundamental equations, 

We know (see e. g. Brown, Bolstreli 1959) that the Wilkinson model of the giant dipole 
resonance (summarized by (2.13)) fails to give the right position of its maximum. Let us see 
what will happen to our R(w) if we switch the interaction (2.10) on. Then it is natural to base 
our analysis on the following equation for K(afyó; w) 


R(affy0; w) =ð. 05, K (aD; 0) +K (xB; o) X; Obur; w)K(uryd; w) (3.1) 


w 
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where K,(@ß; w) is a dressed particle-hole propagator and Q(xfuv; ø) is an effective par- 
ticle-hole interaction which is, in general, Fourier transform of a time dependent inter- 


Ka 


K = E 


Fig. 2 


action. This integral equation is graphically presented in Fig. 2. It contains as particular 

cases all the approximations to the giant resonance problem proposed so far. For instance if 

we sum over all possible *bubbles" (see Fig. 3(b)) of bare particle-hole pairs — a procedure 

which reproduces the most important features of the giant resonance phenomenon (see 

Mottelson 1960) we get the following equation for K (graphically presented in Fig. 3(c)) 
K(afyó; w) = 0,,0,K,(28; w) + Kolap; v) >) V(xB; uv) K(uvyö; w) (3:2) 

uv 
where 


V(«B; uv) = <aulv|vh > 


ER DUNS. 
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Let us keep for a while to this approximation and see how our response formalism works 
in this case. The sum in (3.2) is extended over a finite number of states so we have to solve 
a system of linear equations. In order to have more handy notation we label every pair of 
indices by a number: e. g. a, B==p; y, ð=q; m, ver etc. It is more convenient for our 
purposes not to solve (3.2) directly but to express it through a particle-hole 7-matrix 


Fig. 4 
defined by the equation (presented graphically in Fig. 4): 
T(pq; ©) — V(pq) + X; V(pr)Kolr: 0)T(rq; o») (3.3) 
Then 


K( pq; w) = ð Ko( p; 0) + Ko(p; w)T( pq; w)Ky(q: w) (3.4) 
The solution of (3.3) is: 


A( pq: w) 


T(pg; œw) = Aco) (3.5) 
where 

1—V(py py) Ko( p; 0), — V( py P) Ko( pas 0), .... V(pyq). .... — V(pyp,) Ko(p,; 0) 
thd: t i Vipap) Rok ss 0), 1-7(Pa Pa) Kol Pa: e), 2... V( Pog), > — V(pap,) Kol py; 0) 
— V( PP) Kol ys 0), — V( pa Pa) Ko( Pa: 0), .... V( pq), .... 1= V(p,p,) )و كل‎ Pas 0) 
p-th column (3.5a) 

1-F(P+P)Ko( pis ©), — V pi Pa) Ko( pas €), «+» — V(pyp,) Kol Pn ©) 
KATZ =} (Pap) Kolpi; 0), 1— PaPa) Kol Pas ©), +205 (papa) Kol Pai w) (3.5b) 

— V( PP) Ko( py; 0), — V( p, pa) Ko( Pa: 0), ..., 1— V( p, p,) Kol Pps 0) 


| 

m being the number of pairs of states taking part in the dipole vibrations. Sometimes it 
is useful to present. T as a ratio of two polynomials of the effective coupling constant (kept, 
in our notation, hidden in V(pg)). In this form our solution is completely analogous to the 
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Fredholm solution of an integral equation (see e. g. Tricomi 1957, Courant, Hilbert 1937) 


Aa(P9) + ... + (—1)*-14,(pq)‏ — لوم رك 


T(pg; = 
(pqs ©) He a a. ici 
where 
j yaa) 
A= 2, NE A : (3.6a) 
P1<Pe--- «D, 
V(prPı) ... V(pypy) 
and 
0 $5 PE 
V(P19) V( ps p) Ko(P1: 0) ... V(pypy) Ko( py; 0) 
Apg)- Ð V( paq) V( pa pi) Ko( py: 0) ... V(pa pr) Ko( py; ©) (3.6b) 
fip Xp, 1 : , 4 
V( pra) V(pypy)Ro(py: 0) ... V(pr PR) Ko( py; ©) 


Unfortunately it is rather difficult to give any detailed systematic description of the giant 
resonance based on the above formulas without extensive numerical computations and it 
will be postponed to a later paper. Here we shall merely draw some qualitative conclusions 
from our basic formulas (3.1) (3.4) (3.6). 

The resonant frequencies correspond, as we already mentioned at the beginning of § 2, ' 
to the poles of K(pq; @). From (3.4) and (2.12) we see that we still have Wilkinson’s absorb- 
tion (2.13) introduced into (3.4) by the poles of K, functions. Moreover, the zeros of the 
denominator of T(pq; w) (which satisfy the equation 4(w)=0) correspond to some new 
absorbtion frequencies which do not exist jn a system of noninteracting particles. 

The first (“free”) term in (3.4) gives a contribution to R(w) identical to (2.13). One can 
see, however, that the functions involved in the second "interacting" term of (3.4) produce 
a contribution which, provided our shells are not degenerate, cancels completely the first 
one. Indeed in this case the poles of Ko(p; w) functions are well separated and we can extract 
from (3.5) and (3.5a, b) the *exploding" part of T(pq; c) at, say, p-th pole. From the re- 
lation (obtained by inspection of (3.5a, b) 


A(w) = — K(p: o)A(pp; 0) (3.7) 
we get 


A(pq; w) (3.8) 


1(P4:0) — OAD: 0) 


which shows that T(pq; w) contributes to Q(w) at p-th pole only when g=p. Hence we have 


K,(p: ) T( pp; e)Ko(p; 0) = — Ko(p; 0) 
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and the *free”-term contribution becomes cancelled. This fact does not depend on the 
interaction V(p;ps). We may observe, however, that as soon as we introduce any damping 
into the particle-hole propagator the cancellation is incomplete (see below). In the case 
of degenerate shells the cancellation takes place only if V(p; pa)'s fulfill certain conditions 
(see below). 

Suppose o, is a zero of the denominator of (3.5) (i. e. A(w,)==0). R(w) is given, accord- 
ing to our basic formulas (2.1) and (2.4), by the residuum of S(w) at o=@,. Hence 


E Ko( p; )A( pq; o) Ko(q : of (p (q) 
R(wo) = X - ð(w — w,) (3.9) 
A (Á) 


A 
v v=0, 


where 
f(p) = f(ap) = eh(E, — E„)z(«P) 
Let us note that if 


V(«B; w) = V(Ba; w) = V(«B: vu) = V(Ba: vu) (3.10) 


` 


(satisfied e. g. in the case of separable potentials, see below) we may use in calculations of 
w, and R(w) a simplified form of Ky (defined by (2.12)). We introduce a positive quantity 
ا ر‎ T for EE, 'Then 
E 26 
Ko( p; w) = o? — دقع‎ in (3.11) 
In order to illustrate our formulas let us consider the first (in the sense of Fredholm 
expansion (3.6)) approximation. The excitation energies @, satisfy the equation: 


1= (3.12) 


and (3.9) takes the form: 


Sara بكي‎ 


e = 


Ko) = كف‎ 1 RAWA MIN “a 0(0 — Or) (3.13) 
(w? — £2)? a)” 


As we consistently use the Feynman-Dyson description of particle-hole propagation we 
automatically include in (3.11) and hence in (3.12), (3.13) the diagrams going backward 
in time (see Fig. 3(b)) or, in the language of Brown et al. (Brown, Evans, Thouless 1961), we 
take the ground state correlations into account. (3.12) and (3.13) give an approximate descrip- 
tion of the giant resonance not only for weak interactions but also in the case of arbitrary 
strong but separable potential: 


V(pq) = FD(p) Dq) ˆ (3.14) 
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In this case all 4,=0 (see (3.6a)) except .رك‎ Brown et al. (Brown, Evans, Thouless (1961)) 
give some arguments that one may expect (3.14) to reproduce properly the qualitative fea- 
tures of the dipole vibrations. In Fig. 5 is shown the graphical solution of (3.12) (in the case 
of positive V(pp)). We see from (3.13) and Fig. 5 that @;=@, should be interpreted as our 
dipole collective excitation energy: (R(w,) carries most of the dipole strength of E1 absorb- 
tion and all degrees of freedom contribute coherently to it provided V(pq) are positive?. 


Instead, for all ©(w,) large cancellations take place because of different signs of terms con- 
tributing to the numerator of (3.13)3. 

It may happen that certain V(p,p,)’s vanish and the system of linear equations (3.3) 
may reduce to two (or even more) independent systems. Then the determinant (3.5b) is re- 
ducible and defines more than one collective level. 


2) We accept V(pq)>0 in El absorbtion. For details see Brown, Castillejo, Evans (1961). 
3) One may easily find that R(w,) vanish in the degenerate case (w„= „= e) if the following condition is 
satisfied: 


YP? (p) >, V (pp) = Mf (p)V (pa) (a) 
p y 


24 


which in the case of a separable potential becomes: 


Cp‏ ارم ا رن لا تدر 
b p 5‏ 


and which means that f(p) and D(p) “vectors” are to be parallel. 
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A word concerning the “weakness” of V necessary to make (3.12) work (for non-se- 
parable potentials) is in order here. We know from experiment that w, = 2e. Hence 
VK, is of the order V/e provided we are solving A(w)=0 for .يه‎ Thus the expansion 
in the denominator of (3.6) is in powers of V/e which is very probably a small number (V/e<1) 
in real nuclei. In any case, the use of (3.12) for calculating w, seems to be consistent: it pro- 
duces w,>2e if the order of magnitude relation nVæe is valid (n being the number of 
pairs taking part in the dipole vibrations) which in turn gives V/eæn and justifies our 
first order approximation. We must, however, observe that in calculations of w,<a, (see 
Fig. 5) much smaller V is needed in order to make VK, small. Therefore the use of (3.12) 
is more doubtful in this case. 

One can give an analytical expression for w, in the case of degenerate shells (all e,=e): 


w, =Ver+ F2e X, DXp) (3.15) 
p 


We may briefly summarize these (partly known, see e. g. Brown, Castillejo, Evans 1961) 
results as follows: the free particle excitations disappear completely and a collective absorb- 
tion line, pushed up on the energy scale, carries most of the dipole strength; besides, a host 
of weakly excited lines appear close to the free particle excitation lines. This qualitative 
description seems to be in agreement with experimental results. 


(b) The width. 


Our formulation is very convenient for introducing widths into the giant resonance 
calculations. As was emphasized already we do not start from the first principles and our 
single particle (— hole) propagators (2.11) represent a kind of quasiparticles. Therefore 
we do not make essentially any new physical assumption if we introduce widths for quasi 
-particles and quasi-holes*: 


| 6, (a) 1 0_(a) 
E = (Eg — tye) E — (Ex+ 1%.) 
. This formula, in other words, represents propagation of a “dressed” particle (hole). As I 
want to describe quasi-particles (— holes) in the spirit of Galitski’s paper (Galitski 1958) 
y. y do not depend on Æ and I can perform similar integration as in (2.12) and get: 


G(a, E) (3.16) 


Kalaß; w) = — 0.(8)0- (a) 0 (9)0, (a) 


H= EE + i(yg + Ja) i (ES Eg + Ey — (yet 55) (3.1 


Now we can repeat our previous caleulations (3.3)—(3.6) with Ky>K,. We see imme- 
diately that the A(w)=0 equation has now complex roots and we interpret their imaginary 
parts as the widths of our absorbtion lines. In this way we express them in terms of the sin- 
gle-particle (— hole) widths. As long as we accept (3.1) as an equation describing our 


1 Mes 7 ABAK , A 
) After introduction of widths we do not have a complete cancellation of contributions of the “free” poles 


as geroribed below eq. (3.8). In a real nucleus, however, these resonances will probably merge with those occuring 
at (0,5 
ar 
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dipole absorbtion and use time independent effective interaction (i. e. Q is w-independent) 
Ky > K procedure is the only source of the widths. Then if we know y's we can, in principle, 
calculate from (3.5b) the widths of all excited lines. 

For the sake of illustration let us take again the first Fredholm approximation (or the 
separable potential). Our secular equation is now 


26, — 12 > 
Le > ap y (29) (3.18) 


w? — (e, — ih 
= (£5 — iy, 


where y, is now the sum of single particle and single hole widths and وعدن‎ — il. Taking 
the imaginary parts of both sides of (3.18) we get an approximate equation for I, (c, , €, 


>l,, Yo): 


V( pp) x 


qd X (o? — & a 237 d: 2VpEp 
= (w7 — ep + yp — 17)? + MEpYp — لبه‎ r)? 


20,€, 1 a 
s pr (o8 ety? — T3)? Ay =a)? V (29) | (3.19) 


From this formula we get an approximate expressions for all /"s based on the behaviour 
of e, and w, shown in Fig. 5. For “single particle” excitations 


2 
) | Yptp 
JE R 
1 (oF — eb aE? AE 
b 


*1 , 
€p A Wr 
à E 3 (my 12 A OM Fm) | ŻE: (3.20) 
Instead for the collective width we get 


cx Z V(pp) (e, — ep) + > ep (pp)(w. — &)? 


3.21 
20, X EpV(pp)(@e — £p)? SR 
E 


Ji asus 


where y,, is an average of all single particle-hole y's and we should assume that no y's 
differ very much from each other, this being most likely the case for closed shell nuclei. 
In the degenerate case 

w+ e? 


20.8 


Mom Nor (3.22) 

We should remark that the above formulas are comparatively complicated because 

they take into account diagrams going backward in time. Nevertheless they show that/,'s 

and J’, are only sligthly different from y’s. One may easily check that this is the case for 
very wide variations of quantities which determine Ts. 

As was already pointed out by Brown et al. (Brown, Castillejo, Evans 1961) the 7”s 

for closed shell nuclei should be small because the elementary excitations are from a bound 
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state to a bound state so they have small y’s. Now let us see what will happen if we introduce 
7 7 E= ^c = بح‎ 

a few valence nucleons. Let us take J’, and work it out assuming that all e, we and that there 

are two groups of y’s: a large group of small y’s approximately equal to each other and 


a few large y,’s of the valence particles. From (3.19) we get: 


7, M Vv 7 AF 
V( pp) + — V(p'p | 
2 2 2 (PP) % 3 


waT 


icr 280, A 1 V(pp) + X V(p'p) 
p p' 


(3.23) 


Where I, is the collective width including contributions from the valence nucleons, p runs 
over closed shell nucleons, and p' — over the valence nucleons. To have the crudest estimate 

= ATZ DIEN for n closed s sitions, m valence tran- 
of I, let us suppose that V(pp)==V(p'p'). then for n closed shell transitions, m ve 


sitions 


5 m r 
Degrees (3.24) 
y n 
As y, might be several times bigger than y we see that even a few valence nucleons may 
provide a very essential contributions to the I' of our collective excitations. 


(c) The effective interaction. 


So far we have considered only the exchange vertex (see Fig. 6a) as a kernel generating 
dipole vibrations. But there is another simple possibility of iterating the direct vertex (see 
Fig. 65) and some authors actually included both of them in the giant resonance calcula- 
tions (Brown, Castillejo, Evans 1961). It would seem, however, that although 65 vertex 


y d T d " 


Fig. 6 


may play an important role in the process of absorbtion of the electromagnetic radiation 
by nuclei, it has virtually nothing to do with the collective giant resonance state. 

As the interaction takes place in the presence of the other nucleons, to our time inde- 
pendent interactions 6a and 6b the time dependent ones should be added (e. g.) to 6b — but 
not to 6a — the series of bubbles shown in Fig. 7). Let us consider the graph in Fig. 7. If « 
and م‎ belong to the neighbouring shells the *bubbles" which make them interact with each 


other correspond to lifting a particle two shells up because of the angular momentum and 
parity conservation, 
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The next point is that the introduction of time dependent interactions results in quite 
different modifications in cases 6a and 6b. We shall mark the difference by putting an arrow 
(see Fig. 6): c corresponds to a and d — to b. The introduction of c (which cannot have 
the form shown in Fig. 7 because of the conservation laws) as an effective interaction into 
our equation (3.2) leads us to (3.1) (apart from the “dressing” procedure). Hence nothing 
basic is changed and we have an eigenvalue equation for w. It seems, therefore, in this case 


Erw,a 


Fig. 7 


to be reasonable to use an effective time independent interaction as we did previously. If, 
however, we iterate the kernel 7 we find that there is no equation for K(uByó; w) function. . 


We can only construct an equation for K(aB, E; yó, E’; w) which is related to K as follows 
K(aßyd; w) = | dEdE’K(aß, E; yò, E'; o) 
Namely 
K(ap, E; yd, E'; w) = ö(E — E')6,,05,G,(2, E+ w)G,(B, E) + 


i»! f dE" © بيهام‎ E+ e)G«(B, E)Q(au, Bv; E — E")K(uv, E"; yò, E';w) (3.25) 


uv 
One can of course integrate over E’ but one cannot get rid of E and obtain an eigenvalue 


equation for w. | i 

Without attempting to solve (3.25) we present some indications showing how essentially 
6b vertex may be changed by the time dependent correction of Fig. 7. Let us split our effec- 
tive interaction Q into the time independent V and the time dependent T: 


TEY (3.26) 


We iterate now the kernel (3.26) integrating and summing over all free energy variables 


and «, 6 ... — quantum numbers. Symbolically : 


K= | G,6+ f GGW+ 1) 66, + f CoCo (V+ 1) GGF + 1) GyGy +... (3.27) 
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In (3.27) beside the series generated by V part of the interaction: 
% = | (6,7646, + | G6, 6,6, CoCo + 
we have also series containing /. For instance: 
m= | GoGo! G,G+ [ 6461 606,766+ f 6.166 GG GoGo + ... 
#a = | 6,601 6,691 GG, + f GoGol 1و6,6‎ CoV GoGo + -.. 


etc. Let us suppose that x, contains a certain pole. The point is now that the other x*s con- 
tain the same pole plus some new ones, each of them carrying a certain dipole strength. 
Hence the introduction of a time dependent interaction causes a certain redistribution of 
the dipole strength and tends to blur the single collective state produced by the time inde- 
pendent part of the interaction. 

Let us illustrate the situation for / shown in Fig. 7 in the case of the degenerate model, 
assuming the separable potential (3.14). We limit ourselves to the influence of the x, series 


On Xp. 
<ap|Q(E — E')|Bó> = FD(ayyD(Bo) x 
FAY úc" FASG V | 
i E EL E (3.28) 


Our calculations are for closed shell nuclei, hence the energy dependent (i. e. time depend- 
ent) part of Q is constructed from excitations lifting particles two shells up. Hence A =2e. 
The corresponding collective excitation energy w,>2e since the isospin conservation im- 
plies the isospin triplet oscillation, hence a repulsive potential. 
> = 2 Dim). 
uv 

If we evaluate [G,G„,VG„G„y (vertex 6b) and [G,GILG,G„ (vertex 7) we find that they 

differ by a multiplicative factor 


We 
u . 
w — à, — E+ iy (3.29) 
Thus 
A 
IE Y 2 
M D? (uv) 
xo + x= |1+ Mg = DX, 


© — à, — دع‎ iq 
o, (the frequency generated by V) is smaller than o, +ë. Hence the x, series tends to cancel 
xo at W=W,. If we assume w,=2e we obtain: 


_ 46/0, — 1 
VI 


(3.31) 
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For w 25, 35, 3.32 p= 1, 0.5, 0.1 respectively. We see therefore that the dipole strength 
is partly taken from the pole of x) and shifted to the other poles. 

One can show that x x etc. have also poles at energies 20, +e, 30, +2 etc. Moreover 
one can construct from (3.27) plenty of different series from x, having poles at different 
energies. All this plethora of excitations produces a certain distribution of the dipole strength 
over a wide energy range. It would therefore seem a doubtful procedure to include 6b vertex 
in calcultations of the collective dipole state just because it is so vulnerable and drastically 
changes its role when influenced by the many-body background. 

The author gratefully acknowledges interesting conversations with professor G. Brown 
and professor B. Mottelson on the subject treated here. 
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THE SOLUTION OF THE EQUATION OF ELECTRIC CONDUCTIVITY 
FOR A THIN PLATE IN AN ARBITRARY MAGNETIC FIELD FOR 
A NON-CONSTANT TEMPERATURE DISTRIBUTION 
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The equations of electric and thermal conductivity in a magnetic field are considered for the 
case of a thin plate electrically insulated. The same temperature distribution is applied to the 
upper and lower surfaces of the plate. In this case the set of equations may be reduced to a single 
equation of only two variables, as we may admit that the temperature is constant throughout the 
entire thickness of the conductor. Thus, it is not necessary to solve the equation of thermal con- 
ductivity. A vector potential is introduced and an equation for it is derived. This equation is 
equivalent to one mentioned above, but is more easily solved. The special case of a rectangular 
plate is discussed. It is shown that the stationary electric current gives rise to a deviation from , 
the Nernst-effect in an insulated conductor. 


1. Introduction 


In a previous paper by this author (Stachowiak 1961, cf. also Stachowiak 1960a, b) 
equations for the thermal and electric conductivity of conductors in an external magnetic 
field have been derived in the linear approximation, in terms of the latter. These form a set 
of two partial differential equations of the elliptic type. This set has been discussed with 
respect to the boundary conditions corresponding to the electric insulation of the conductor 
and some prescribed distribution of the temperature at its surface. It has been proved that 
the solution of this problem is unique, and a general method of finding this solution has 
been given. 

The purpose of the present paper is to give a discussion of the mentioned equations 
for the case of a conductor having the form of a thin plate. To the upper and lower surfaces 
of this plate we apply identical temperature fields. If the plate is sufficiently thin, we may 
assume that the temperature does not change across the plate and has the same value as at 
the surface. Therefore, it needs not be computed from the set of equations. We reduce the 
problem to solving only one partial differential equation in two variables for the electro- 
chemical potential with a boundary condition for the normal derivative along the side border 
of the thin plate. Making use of the fact that the divergence of the electric current vanishes, 
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we introduce a new potential whose curl is equal to the current. The equation of this new 
potential is simpler than that of the electrochemical potential. After the general equations 


have been derived some special cases are discussed. 


2. The equation of electrochemical potential for a thin plate 


In the presence of a magnetic field, and for an isotropic conductor, the electric current 
density may be expressed by the formula (cf. Landau & Lifshits 1957) 


= J= LVu+ ha” T+ kV u, B] + Lu[VT, B], (1) 


wherein y is the electrochemical potential, B — the vector of magnetic induction, T — the 
temperature, and ly, la, lg» ha — Onsager's coefficients. In formula (1), terms non-linear 
in B have been neglected. We assume that Onsager’s coefficients depend only on the tem- 
perature. 

In a previous paper (Stachowiak 1961), we have proved that, if the normal compo- 
nent of the electric current vanishes at the surface of the conductor (i. e. the conductor is 
electrically insulated), yu may be written in the form 


PE ŚR bróð (2) 


where Ap is a quantity proportional to B. 


Substituting (2) in (1) and neglecting the terms nonlinear in B, we obtain 


— J=1V—L{VT, B] (3) 
wherein 
_ dhs — hha 
L= ee: (4) 
Let us consider a thin plate situated in the plane XY, and assume that at its surface the 
conditions 
J, = 0, (5a) 
I(x, y, 0) = T(x, y, d) (5b) 


are satisfied. Here J, is the normal component of the electric current, and d — the thickness 
of the plate. The condition of stationarity is of the form 


V-J=0. (6) 
We now average equation (6) over the thickness of the plate, making use of the condition 
(5a): i : 
E V. Jdz = NIĘ +2) dz = 0. (7) 
0 0 
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We substitute Eq. (3) in the foregoing equation and take into account the condition of Eq. 
(5b). Applying the mean-value theorem, we obtain 


(er on (2222, seen) 
9x? dy? dx Ox dy dy 
(2 oT | BR. z 

dx dy dy Ox] ا‎ 0 © 
assuming that the magnetic field is constant throughout the entire thickness of the plate. 4 
is the derivative of J, with respect to the temperature. The functions on the left-hand side 
of (8) are taken at a point zy(x, y) within the plate. If the latter is sufficiently thin, the quan- 
tities on the left-hand side of (8) do not differ considerably from those at the surface. There- 
fore, we shall henceforth understand the functions in (8) to be taken at the surface. This 
enables us to avoid considerating the equation of thermal conductivity. Substituting (3) 
in (5a), we obtain the boundary condition 


A 
ŚW اا‎ at the edge of the plate. (9) 
On h ot 
3 BE Ð 1 NA ee tod 
Here, T denotes the normal derivative and nir the tangential derivative with respect 
n 


to the border. 
Eq. (8) and condition (9) define a boundary problem of Neumann’s type for Ó. 


3. Introduction of a vector potential for the current 


As the problem is two-dimensional, the curl of the current has only one non-vanishing 


component 
= OJ, 9J, 10 
up eee Q0) 


Here, J, and J, are, of course, the components of the current averaged over z. Substituting 
(3) in (9), we obtain 


1 387 (637 
Iv. Tie e - 1, (55 | za) | 
p [87 2p _ 9 
wire Ror 


+€ 


7 9B, oT | 8B, = 
"NOx 0x dy dy]? 
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where L/ denotes the derivative of L with respect to the temperature. Using (3) and neglecting 
terms quadratic in B, we eliminate derivatives of @ from the foregoing equation and obtain 


2 2 s IT 27 
V A le | 


yy 


9 2 er 
eJ. = 44 = 0 (7 ) 
Ox dy 


= RE (12a) 
e 

^ 

A 
0 (12b) 


aT)? aT)? OB, OT ƏB, OT 
- +L = + z 13 
3 IÐ | EJ | ne E: x | oy 2 (13) 


= 0 along the edge of the plate. 


The arbitrariness of the additive constant in A resulting from the definition of eq. (12) allows 


us to assume that 
A= 0 at the edge. (14) 


Eqs. (13) and (14) define a boundary problem of Dirichlet’s type for A. 


4. Special cases 


In Sections 2 and 3 we have obtained two equivalent boundary problems for the poten- 
tial and the electric current in a thin plate. Solution of the problem (8), (9) yields the distri- 
bution of the potential, whereas (13) and (14) is a problem leading to the electric current 
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distribution. The relationship (3), of course, admits of the immediate transition from the 
one problem to the other. With regard to the convenience of either method of approach, 
the problem for A is mathematically simpler, as only the equation is non-homogeneous 
(the boundary condition (14) is homogeneous). If the temperature dependence of Onsager's 
coefficients is neglected, Eqs. (8) and (13) assume the form 


dwa op L (eB, oT AB, OT (8) 
Ox? | شيرق‎ L \ dx dy dy Ox]' 
034 3A SM CUT. OB, OT OB, dT 
EX: pe ل‎ te | | 5 | 7 13 
9x? 9y? ur = | x (2 Ox oy Oy ) ال‎ 
If the component B, is constant, we have 
Bp 9?g 42 
— =Ü 8 
943 * Oy? ©) 
A PA JR a Jed 7 
SV ل ا‎ (eH | = Mae 


We shall consider this special problem in greater detail in order to demonstrate some 
physical facts resulting from the equation of electric conductivity in a magnetic field. If 
Onsager's coefficients and B, are constant, we see from (8") and (9) that q depends only 
on the values of the temperature in the edge of the plate. 


Let us consider the rectangular plate shown in fig. 1. For the sake of simplicity, we 
assume the following temperature distribution along the sides of the rectangle: 
T(x, 0) = T(x, a) = T(x), (15) 
JUL cit A 
T(, y) — T(). (16) 


By the foregoing boundary conditions, the problem of Eqs. (8) and ©) can be solved di- 
rectly (cf. Morse & Feshbach 1953 a); in order, however, to avoid a Fourier expansion which 
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would complicate and render less clear the essence of the problem, it is more advisable to 
consider the equivalent problem (13°), (14). In this case the solution may be written imme- 
diately, as the Green function of this problem is well known (cf. Morse & Feshbach 1953 b): 


3 PUT COMI A 
A(x, y) =- LX - oP Dye 23) G(x, y|xo, Vo) Exod yo, (17) 


where S is the surface of the rectangle, and 


TVS E EET 
B ا ل الل‎ 
2 a a 
G(x, y|Xg, Yo) = ) 
: 1 zn ¢ 
naa n sh — 21 
a 
zn(l + x) an(l—x) - 
sh sh 9" for % > x, 
j a a 


E an(l + xo) | an(l— x) Re 


a a 


(18) 


In this case, however, it is simpler to proceed as indicated in the Appendix. 

Let us note that the function y depends only on the temperature along the border. If 
we compute p by means of (3) and (12) from the function being the solution of Eq. (13) 
at the prescribed temperature distribution (15) and (16) along the border, we may assume 
that T is independent of y. This assumption may have some effect on 4 and on the current, 
but not on ©. Knowing g and the actual temperature distribution within the rectangle, we 
can find the actual current by means of (3). Thus, the solution of Eq. (13") is easily obtained 
from the solution of 


E Le, )19( 


assuming that the temperature distribution satisfies the conditions (15), (16). 
By results derived in the Appendix, easy computations yield 
a(2m +1) 
n — — 


si y 
Aa a : 
A(x, y) 2 —LB, - 
(% y) s, X az? (2m +1)? | (2m + 1) 9l ii 
sh — 2 


a 
SALA 1 +1) aT | x(2m-+1) 
x |»? ca fa " (1 xo) dx, + 
HE dT | 2m +1 
NUM US EET. 


259 


We now compute the voltage between the corners of the rectangle: 


u(l, a) —u(L, 0) = p(l, a) —g(i, 0) = Ma dy. (21) 
A كرالك‎ 
0 
By the relationship of Eq. (3) and the definition of Eq. (12), 
op Jy  LB,dT ١ LB,dT , 194 
dy 2 IUE dx DNS TINO 22) 


Substituting (22) in Eq. (21), we have 


Bain 1 94 
l, a) — u(L 0 = i 7 
>) |= DES da 4 ki fÐ) o. (29) 
0 


The second term on the right-hand side of the last equation represents a deviation from the 
Nernst-effect caused by the existence of an electric current. We compute it on the assump- 
tion of the following temperature distribution: 

T(x) = Ty = e|x|, (24) 
where c is a positive constant. Physically, such a distribution is not strictly possible, but it 
simplifies our computations. From (20), according to this assumption, 4 is seen to be little 
different from the values it assumes in physically possible temperature fields, i. e. such that 


the temperature maximum is less sharp. 
By the distribution of Eq. (24), we have 


dz | 
dat = — 2có(x), 
and, for x >0, 
a geam 
A(x, y) = 4cLB,a » mi D BTE sin A y. (25) 
0 ch — 
a 
Substituting (25) in (23), we obtain 
LB = 1 
L, a) — u(l, 0) = = 26 
N ~ GT D p am), 0) 
= h si 


If 7 is not much smaller than a, the first term of the series in. brackets yields the deviation 


to a sufficient degree of accuracy. 


5. Consideration of the temperature dependence of Onsager's coefficients 


In the general case, neither in Eq. (8) nor in (13) can the variables be separated. In order 
to make this separation possible, we must choose a temperature distribution depending 
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on one variable only. Let us assume the distribution 
لات ل‎ (27) 


Eqs. (8) and (13) can now be written as follows: 


og. a 140,0 L3B, dT 
Bp 3p | l dT dg + dT A 
Qu Gaye P da) Li S 
U LR WAW NIA 
ga. det rr e 
mo dt\* 3B, dT E 
TEDA LLL (1377) 


In this case also the problem for A is simpler than that for g. Therefore, only Eq. (13) 
will be considered. The problem (137), (14) is particularly simple if we assume B, in the 
form 


B, = f(x) g (y). (28) 


when we may use the method given in the Appendix. However, on separation of the varia- 
bles, we obtain for X(x) a differential equation with variable coefficients for which the so- 
lution in general form cannot be written. To avoid this difficulty, we assume the temperature 
distribution (24). Moreover, we assume that cl is small and Onsager’s coefficients depend 
weakly and linearly on the temperature, i. e.: 


cll, €l,, cli'«L. (29) 


Hence, we can write 


L(T) = HO + (T — TONG. (30) 
L(T) — L(T) + (T — ToL (Ty). 
In particular, by the distribution (24), we have 
L(x) = 4) — c|xll(T„), 
(31) 


L(x) = L(0) — c|x|L'(T)). 
On taking into account the assumptions made above, Eq. (13%) assumes the form: 


A 2A A(T) ? 94 c 
du + yr * hn) CBM mo = gl) | - 2619/6) () + 


"du h (T, j 
+ c? L (Ty) — Lus un) | f(x) — esign x L(T)) f' (x) + cx Lo. (32) 


The solution of Eq. (32) for the boundary condition (14) can be represented by the integral 


1 
Alx, y) = f hæ) T's, ys xdg, (33) 
-1 


761 


wherein J” fulfills the equation 


O PA rae he oe ar 
2,3 = ae LT) CHER يد‎ en g(y)ð(x — x) (34) 


and vanishes along the sides of the rectangle. 


h(a) = UT) — e Ka ue LT) | 9 + 


T LT ye sign xf '(x) — c*?xL'(Tq) f' (a). (35) 


is readily computed by means of the results obtained in the Appendix. 

If B. (x, y) is not of the form (28), we have to compute the Green-function, which in the 
present case, presents no great difficulties. 

The author wishes to thank Professor R. S. Ingarden, Professor J. Łopuszański and 
mgr W. Zietek for their valuable advice and helpful discussions. 


APPENDIX 
Let us consider the differential equation 
Û (x,y) = — FE) (1) 
defined in the rectangle of Fig. 1, with the boundary condition 
p= 0 along the sides of the rectangle. (IT) 


Óis a second order differential operator of the elliptic type in two variables. The solution 


of the problems (I), (II) is of the form 


y», y) = h G(x, y | xo. Yo) F (xo) (V0) dodo » (HT) 
5 


wherein G(x, y|x, Yo) is the Green function of the problem. 
Moreover, let us consider the analogous equation 


Ô T(x, y; x’) = — &()o(x — x’). (IV) 


Its solution is of the form 
D(x, y; x) = f dygg(yo) G(x, Fir» yo). (V) 
0 
This enables us to write, on substitution of (V) in (HI), 
yx, y) = j/ dxo f (x9) T (x, ¥; xo) - (VI) 
e] 


T(x, y; xo) is the potential of a single layer, as may be seen from Eq. (IV). This potential is 
easily obtained, provided the operator Ó allows to separate the variables. Let us assume this 
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to be the case. The straight line x=x, divides the rectangle into two smaller rectangles in 
which the following homogeneous equation holds: 

OT (x, p; %) = 0. 
On separation of the variables, we obtain J'in the form of the series 


I (x, yix) = DTR (x) Xg (a) Y„(y) for x > 29. 
$ (VII) 
T(x, y; xo) = Le G9) Xe (x) Y4(y) for x > x. 
k 
Herewith, the following boundary conditions are fulfilled: 


AD = X; (-1) 2 Y,(0 = Y,(a) =0. (VIII) 


We obtain the coefficients IF (xy) and Ty (x) from the following conditions on the single 
layer: 
1) continuity of the function 


P(x FO, 7: 3%) = Ix — 0,y5 29); (IX) 


2) discontinuities of its derivative (cf. Miranda 1955) 


9 AJD 
= - go). (X) 
CX |x,+0 CX ix —0 
We now expand g(y) in the series 
gly) = Dave YY) - (AD 
k 


On substituting (XI) in (X), we obtain from (IX) and (X) a set of two equations for TĘGI 
and f(x) for all k's. Hence, 
Ds zak (o) | 
, Xp (20) Xe (Ag) EAR (29) Xk m) 
E Xe (x f 
I (a) ee ال‎ (XI) 
Xx (xo) Xx. (x9) — Xk (%9) Xk (2) 


The function y can be obtained from the formula 


* L 
EIER TERN EEE TER Rd], X 
-l x 


which was obtained by substitution of (XII) in (VII) and of (VII) in (VD. 


lf the right-hand side of the equation for y is not of the form f(x)g(y), we can obtain 
Green's function by employing a method very similar to the above, presented in the book 


of Morse and Feshbach (1953 b). 
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PROJECTIONS AS OBSERVABLES 


By JERZY PLEBANSKI AND JOANNA RYTEN 


University of Warsaw, Warsaw, Poland, and Institute 
of Physies of Polish Academy of Sciences 


(Received April 12, 1961) 


A pair of events consisting in the emission of lisht pulses may be observed by observers 
located at a plane at infinity (determined by the orthogonal direction to it). The data concerning 
our events obtained by these observers can be considered as observable quantities independent 
of the coordinates used in the fmite regions of Riemannian space-time. This general idea (first 
given in [1]) is discussed and supported by explicit computations in the case when the linear 
approximation applies. In particular, the general formulae for projections of a “distance” between 
two events and for their projected “temporal separation” are derived. With the help of these 
formulae the projection of the geodesic motion is studied. Eventually the formulae obtained are 
applied in the case of a field induced by a rotating body. 


Introduction 


It is well known that the theory of measurements in general relativity does not almost 
exist. We are able sometimes to derive mathematical answers to which the theory leads; 
unfortunately very often we do not know what is the physical content of them. Indeed, so 
long as the results derived depend on particular coordinates, there always remains the ne- 
cessity of finding out whether these results are objective, i. e. whether they may lead to ob- 
servable effects. 

There have been several attempts to construct a theory of observables in general rela- 
tivity. One can try to attack the problem according to one of the two following obvious ap- 
proaches. First, we may search for some objective meaning which can be attributed to the coor- 
dinates in which we have solved a given physical problem. In other words, we may try to 
establish the physical and geometrical meaning of that system of reference to which the 
coordinates x“. refer. This program may be carried out in several special cases when, because 
of some additional circumstances (e. g. some invariancy under some group of motions), 
we are able to give definite prescriptions how to measure these coordinates. 


1 The harmonic coordinates so convenient and elegant in mathematical applications, unfortunately suffer 
an essential disadvantage: namely, even the most enthusiastic supporters of these coordinates (as the only 


objective ones) were unable to give a prescription how to measure them. 
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The second possible line of attack consists in dealing with directly measurable qa 
which, from their very nature, cannot depend on the coordinates used in computing their 
values. 

More formal approaches, are possible also like that of Pirani [2], who proposes to con- 
sider as physical quantities the projections of tensors on some tetrad frelds. 

The paper of Newman and Goldberg [3] which methodologically corresponds to the 
second possible approach discussed above forms, no doubt, some step forward. The authors 
discuss there the mechanism according to which we may construct in general relativity such 
observable? quantities as luminosities, angles etc. Altogether, however, the existing situation 
is still pretty confused; by no means we could say that we possess already a satisfactory theory 
of measurement in general relativity. 

The difficulties in question are especially crucial from the point of view of the theory 
of the equations of motion. This theory is able to derive post-newtonian equations, i. e. 
newtonian equations corrected by terms proportional to c*. These corrections, however, 
do depend on the coordinates used. Having post-newtonian equations and an approximate 
metric in arbitrary coordinates, we would like to know to what extent these results correspond 
to reality and to what extent they reflect spurious effects caused only by an unappropriate 
choice of the coordinates. 

In the work of Infeld and one of us [1] in connexion with this problem a simple method 
was proposed: this method consists essentially in the simple remark that these things which 
an observer, in particular an observer at infinity, sees (in the literal sense) being automati- 
cally objective cannot depend on the coordinates used in the computations. But observers 
at infinity have a great advantage relatively to those at finite distances; communicating be- 
tween themselves they can use concepts of the pseudo-euclidean geometry of flat space-time. 
That all leads straightforwardly to the following idea: to every given motion in Riemannian 
space-time and to any direction determining a plane at infinity orthogonal to it we may 
attribute what we may call the “projected motion” or the “shadow motion”. 

Consider a time-like world line in a finite region of space-time x*(s); now let /7(n*) be 
the plane at infinity orthogonal to the cartesian vector n“. (We assume that our world is asymp- 
totically flat). Let the point which moves along a*(s) be sending permanently light rays in 
all directions. Now, let us select out of all light rays which were emitted from a*(s) in the 
“moment” x(s) that one which will arrive at infinity with the direction n* (of course after an 
infinite time). Let us repeat this construction starting this time from the event a*(s) + dx%(s). 
Obviously (see Fig. 1) the distance between the two points on 77 in which our rays reach it, 
dX „is well defined; so is the interval of time on *Minkowskian" clocks of the observers at 
infinity consisting in the temporal separation between the arrivals of both our rays on II. 
(It takes an infinite time for our rays to reach 77 but the difference of these two infinities 


* In the literature concerning the problem of the quantization of gravity the expression “physical obser- 
vables" is often used, see e. g. [4], [5], [6]. These “observables”, however, have only very little in common with 
quantities which can be directly observed and henceforth are independent from any coordinate system in which 
we are leading computations. These quantum observables correspond merely to a set of independent — from 
the point of view of the canonical formalism — dynamical degrees of freedom; the objective character of their 
definition from the point of view of the covariancy principle seems to be doubtful. 
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determines uniquely a finite quantity dX2). Therefore it is clear that if we have a succession 
of events at finite distances forming a world line, one can attribute to it with the help of our 
construction a “projection”, a “shadow motion” on every plane /7 consisting in the depend- 
ence X1. (X9. n). These quantities are obviously independent of the coordinates in terms 
of which we labelled the events of our world line and of the light rays. If the curvature tends 
to naught at infinity, and if the allowed coordinate transformations preserve the cartesian 


Fig. 1 


character of the coordinates at infinity, we may say that under coordinate transformations 
our light rays behave like rubber lines with ends rigidly fixed to the plane at infinity. Coordi- 
nate transformations may change the shape of these lines; but, from the point of view 
of the observers at /7 only their ends matter. Thus, so to speak, the relationships between 
these ends form an objective physical information about the moving bodies in finite parts 
of Riemannian space-time to which our *rubber-like" light rays are attached. 

Hence, there is no doubt that the set of shadow motions for all directions n“ forms ' 
a set of objective informations about the motion in question. Moreover, in the case of a flat 
world these informations are sufficient to determine completely the motion itself. We are 
not claiming that these informations are all objective informations which one can derive 
about the motion. Even the explicit knowledge of all our *shadow motions" could not satisfy 
us; it would be much better to be able to avail oneself of observables concerning these 
things which observers located at finite distances could perceive. 

Nevertheless it seems reasonable, even if we had actually at our disposal only our *shad- 
ow observables”, to develop a little deeper their mathematical theory and to inve- 
stigate its consequences. This is just the aim of the present paper. 

In [1] the projection of the motions of a test particle in the field of a heavy particle as 
well as the projection of a motion of two heavy particles were performed. This was done, 
however, in the case of solutions obtained by means of E. I. H. approximation method and 
all the formulae were derived only in the case when this method applies. 

In the paper of one of us [7] there were given formulae for the world lines of light rays. 
The only assumption made about the gravitational field was that the deviations of the metric 
Ag” =g” — n” (gł = (— g)” g”) from Galilean values 7” are small and vanish at infinity 
together with their derivatives at least as O(r ^), r = (x*x*). The formulae for the light rays 
were given with an accuracy up to O(4g*). In this paper we shall apply these formulae in 
order to compute the projections also with an accuracy up to O(Ag^); the results in this appro- 
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ximation seem applicable in most of the physical problems we are dealing with. The question 
of finding analogical formulae with better accuracy (e. g. up to O(A?g"")) would be a simple 
generalization of the presented procedure. 

In section 2 we give general formulae for the distance X4 (x^, °, n?) and for dX (x^, 
a9, n“). The procedure of projection is as follows: first we take a plane orthogonal? to n“ at 
a finite distance from the point x“, and we project it on this plane by using a light ray. Then 
going with the “arithmetical” distance between x“ and the plane to infinity, we derive the 
projection on the plane // at infinity. 

In the last part, we use our method to project a geodesic motion in the *weak" gravita- 
tional field; by “weak” field we understand a field corresponding to assumptions under which 
the world lines of light rays with the accuracy up to O(Ag“’) were derived in [7]: geodesic 
motions are here considered with the same accuracy. 

As an example, we discuss finally the case of an gravitational field induced by a slowly 
rotating body. 


2. General projection formulae 


Throughout this paper we assume that 
"نو‎ = (— gig" = n" + Ag" (2.1) 


where the quantities Ag” are small in comparison with 1. In [7] it was shown that with the 
accuracy up to O(/9/) the world lines of light rays leaving the point x“ at the time 
x? have the form? 


X? 
XX) = a^ K[X? — 29 — 1 [ dx Ag" (°, of + M(x’? — 2%) keke] 
+ (6% — kk?) [4 | IX — x'9A g^r(x'9, x* k (s — x9))k-k- 
J , : v m 
x? 
+ | da'%A g” (x'9, ać + d(x" — Jk] (2.2) 


where k* = [1, k*] and k* is an arbitrary unit vector. (Indices with bars are indices lowered by 
n metric, e. g. k; —n, k^). 

From all our light rays (2.2) we must choose now this one which, while reaching the 
plane orthogonal to a vector n being at a distance / from x“, has the direction n“. If we 
denote by X? the time needed for the ray to reach our plane, we get the connection 
between n* and k* (for detailed calculations see [7]) 


if == n* — (6% =n*n') [4 GSO: x + 0X? — x9)n- 
Xi 
+4 l dx ^A goa, af + n*(a'0 — 2))n-n-] (2.3) 


3 7? e 0 0 0 0 
“Orthogonal” in the usual sense: x@ are treated as simple cartesian coordinates; this orthogonality makes 
covariant sense only in the limit when the plane tends to infinity. 
4 This corresponds to equation (4.11) in [7] in a different notation. 


769 


so the light ray which after the time X? has the direction n^ is:. 


لوم 
AX"(X9) 7 x^ + n* + [X0— x9— 1 | dx A g(x", x* + n*(x'9 — x9)) nn | +‏ 


x? 


+ (8° — nn?) (x? — X°) [A g^ (XP, + (X — 29) n; + 


XI 
+4 " dx A gy (x'9, x + n (x0 — x°))n-n.] xm 
x? 


ze 
+} f dx — 4) Agr (x0, + n? — æ))n m + 


X9 
1 if dx! Ag” (a6, s + n (0 — 39) n. (2.4) 


it has therefore the form 
AUX?) — x n*X, + X] (2.5) 


where the explicit formulae for X and X may be immediately written down by comparison 
of (2.5) with (2.4). Now, from the condition that Xy(X)=1 we can find X? 


X? = I+ 40+ 122 r DE x?))nn- : (2.6) 
E. 
Differentiating at constant | we get dX? 
dX) = dx? + 4 [A g” (x9 + L x* + n'l) — Ag" (x°, x) |n-n-dx® + 
xl 
+3 ji ONE r (2x0) nn (do? — n dx). (22) 
» 
Inserting (2.6) into the expression for X^ (X?) we get 
a (XP) =. (ô? — nên?) [— LA g” (I+ x9, x + n'ln- + 
Kuil 
+4 T Date VAN A en للا‎ w0))n.n. + 
0 x+1 


x il dx'* A g” (xo, x! + n) — a?))n-] (2.8) 


and after some simple computations 


SET 
X4 (XP) = (6? — nn?) J dx [A gts, E 0 — x°)) 
x? 


— &A gt? (x, af + n(x? — a9) Jnn (2.9) 
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= e ee : TS ; 
This together with (2.7) gives us the projections of our event (x^, x?) on a plane at a ASRR l 
from it and orthogonal to n“. If 1 > co this plane becomes the plane II at infinity. Let us 
put lim X) = Ks la X^ (X?) = X4 .. After the limiting process, taking into account 


sód 1+00 
wał g"(x0, 45) >0 (uniformly from the point of view of x") for (a)? >œ and 
integrating (2.7) by parts we get our final formulae of projection: 
oo 
dX, = da + 1 dx“ | d'?Agn(^6, x + n'(x'0 — a?)nn- (2.10) 
t 


Xå =4(8%— nn) f ° — ')[24 g(x, af + na — x9) — 


m g^ ^ (7o, a + nî (a^ — a9))n-a-] : (2.11) 


These formulae may be used to establish the correspondence between an arbitrary motion 
in the finite region and its “shadow” at infinity. 

It seems to be of some interest to draw attention to a general conclusion which can be 
drawn from it. Suppose there exists a coordinate system such that metric g“ is stationary 
in it. In this case (2.10) may be written in the form 


oo 
dX& = da + 1 dać | dx'%A gr; (x? + n'a!) ncn. 
0 


oo 
= d[x°+ 3 li dx'°A g*(ax + n*x' 9) n.n] : (2.12) 
0 

From this it is evident that, in spite of the fact that the time dependence for the “shadow” 
is not the same as that for the projected motion, the projection of every periodic motion of 
period Ax? is also periodic and its period AX9 = Ax°. Thus the periods Ax? of the motions 
in finite regions are objective. Every such motion, however, may be regarded as a clock | 
measuring the time x°. Thus this time is also objective: its rate of flow in every finite point 
may be perceived by observers at infinity. 


3. Projection of a geodesic motion 
Let now a*(s) along which our point is moving be a geodesic line, i. e., 


d?x« a do do 
us + Tb "W > 0. (3.1) 


We decompose the initial position x^(s) into its parallel and orthogonal parts with respect 
to n^: 


a= nn) +24, | nx =0. (3.2) 
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Having the differential equation of the geodesic (3.1) and our projection formulae (2.10), 
(2.11) and assuming that in all computations we are interested only in the results with an 
„ accuracy up to O(A g”), we want now to derive explicitly the law which governs the shadow 
motion, namely a law of the form 


dR NIS, 
dX = F = dA: nd (9:3) 


where R] = X] „ + x is treated as a function of X9. and EU X? are related to x*(s), 
xs) accordingly with our projection procedure. 

Intending to derive such an equation, first of all, one should act on (3.1) with n“ or with 
(6% — nn”) and by decomposing x^ accordingly to (3.2) obtain a pair of equations for 
d*x, Px 
D ds? 
respect to ds by those with respect to dx? and subsequently, using (2.10), by differentiation 


respectively. In the next step, one should replace differentiations with 


with respect to dX. In all these steps quantities of the order O(4?g/") should be con- 
sistently neglected; thus, in the internal arguments of 4g”, as well as in the integrants 
and lower limits of the integrals in question, x4, may be replaced by Rj and x? may be 


der, 
identified with X9. Now, the second derivative us may be computed as the sum 
oo 
Xi d?x^ 
of —— 3 and 723. The first term may be derived straightforwardly from (2.10) and 
GA x dA 
(2.11), the second has been already found. But these quantities still do depend on xj; we | 
dx 
know, however, also the second derivative of this quantity xv in terms of X%,, R4, x 
and FF Eventually, therefore, we have to do with a simultaneous system of two 
oo 


differential equations of the second order for the quantities R4 (XL), x (XS) from 
which, however, only the first one has a direet physical meaning. If by integrating the 
second of our equations the non-physical xy (XS) might be eliminated, we would obtain 
the law of the shadow motion in the desired form, 

All the steps discussed above were explicitly performed. The differential equations 
derived are, however, in the general case rather lengthy and not very instructive; we decided 
therefore not to give them here explicitly and to restrict ourselves only to an explanation of 


the whole procedure. 


4. Projections in the case of the rotating body 


As an example let us consider the gravitational field induced by a rotating body. We 
assume that the speed of rotation is small and that the body has spherical symmetry, the 
deformation caused by rotation being negligible. We also assume that the density of our 
rotating sphere of radius / is constant and that the body rotates uniformly with the angular 
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: A sað سر‎ re, 
velocity .كن‎ It was shown in [7] that under these assumptions Ag” is: 


4km 1 4km [3 1 l |e|? l 8911 
ME E Dee 
2k x? 2k res Sy l 3 dul il P 
A 508 > 3 E abc le? Je E € abc x | 9 B 9 5 ar e(l lae) 
A g^ —0 (4.1) 


1 
where m = an (g = const) is the total mass, كل‎ = [escent G s Coh 


is the angular momentum and O(u) is 1 for u > 0 and 0 for u — 0. 
After some simple computations, we get from (2.10) and (2.11) 


2km < 2km l 2k 1 
AU -0 -a „sh gar EA. — — — + ac 
re | © bia) dr a & ; r(r + x) 
> ok n^ Jex“ يو‎ 1 [2 "A 1 xj = 
e Enbc AH eed Á aA 5 r 3p 
9]. 
» = Enon in" Jext, H + dx? A*(3*, n*)O(I — r), SA 
© r 
p 2km x 4k 1% 2k anb Jes =m 
Are Fl = + x Æ Cab. JA ik " ZĘ p Esc Jex L r(r + xj 
27 ai 2k cd 1 
z Æ Gakx1J r(r + xy) Fü Æ E abc J = = r +q] m 
axe — 
+ W Gaj, LEZA, peed, ر‎ (4.3) 


rx), 


where x“ is decomposed according to (3.2) and r?= a%x? = uy T aj. 4° and B“ are known but 
rather complicated functions which can be neglected if we restrict ourselves, as we shall 
do in what follows, to the light rays which do not penetrate through the rotating body. 


In order to obtain the physical interpretation of these formulae, we introduce two unit 
a 


x O 
vectors: u = + and w = Cau n^. Together with n? they form a “dreibein”. 
XL 
Let J? = gn? + Bw“ + yu“. Thus (4.2) and (4.3) can be rewritten as follows: 


. 3 89,2 
dX9,— dx 4 dza Ë Ae SEL 2EB „7 + af 2a 
c r(r+ xj) 0 xn 
2km „1 , 2kB xj , 2ky Mane 
LIPSII w AE الس‎ CH 
2 2km % 6م21‎ 1 2ka x 2k 
X lx = r BEN API ع‎ baise in sb = 6 1 
d m ajj e Soo ad Xj 03 r(r + x) 4 r(r+ xy) (42) 
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The terms in (4.5) and (4.6) which are proportional to m give the projection formulae in the 
Schwarzschild case. The existence of the angular momentum (or ratber of its part ß in the 
direction w”) introduces some corrections into these formulae. In (4.5) there is also a term 
proportional to the y part of J^ but the part of J“ parallel to n does not appear. In (4.6) in 
addition to the component of X4 „parallel to u* we get another one, orthogonal to it and 
parallel to 20%; all the three components of J“ are present. 

The authors would like to thank Professor L. Infeld for his kind interest in this work. 
We are indebted to Mr. B. Mielnik for many illuminating discussions. 
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NUCLEAR ELECTRIC QUADRUPOLE MOMENTS OF 12Sh AND 1223Sb 
BY HUES 
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The hyperfıne structure of the line 4 3030 A (5 p? 2Py,—5p*6s *Ps),) Sb I was investigated by means of a Schü- 
ler hollow cathodedis charge tube as light source and a Fabry-Perot interferometer with multilayer dielectric 
mirrors as high resolving power instrument. Hence, the quadrupole coupling constant was determined for the 
5p26s AP level of the Ys nuclide and, with the electric field gradient as evaluated in the present paper 
for the configuration ps, the electric quadrupole moment of the 124Sb nucleus was obtained. The quadrupole 
moment of the 1?35b nucleus was computed from the ratio (2,93/Q,91, which is known to a high degree of accura- 


cy. The results, without Sternheimer's corrections, are: 


Q(121 Sb) = (— 0.54 + 0.08) x 10-24 cm? 


Q(123 Sb) = (— 0.69 + 0.10) x 10-24 cm?. 


Introduction 


Investigation of the electric quadrupole moments of atomic nuclei is a source of valuable 
information in a number of problems relating to nuclear structure, in particular with respect 
to the exact assessment of the ground states of nuclides. This is essentially the case of the 
nuclides in the neighbourhood of the nucleonie magic numbers. Antimony, which possesses 
one proton outside the closed proton shell (Z —50), is an element of this kind. The electric 
quadrupole moments of either stabile isotope of antimony *3;Sb (I=5/2) and F (LE TÐ) 
have been investigated both. theoretically and experimentally by a number of authors during 
the past few years.! 

It should be stated, however, that the values of the quadrupole moments of either iso- 
tope, as obtained by the various authors, diverge considerably. This is equally true of the 
results obtained by hyperfine structure (hfs) analysis of the spectral lines of Sb and of those 
derived. from microwave measurements in molecules containing Sb nuclei. From the data 
hitherto published, it is hardly possible to explain such considerable divergences, all the 
more so as the hyperfine structure data were obtained by various authors who analyzed the 
energy levels of the same Sb II electron configurations. 


1 A list of the results published hitherto is given in the final part of the paper. 
(775) 
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For this reason, the present author undertook to carry out new measurements of the 
quadrupole coupling constant with regard to the hfs analysis of lines of the Sb arc spectrum, 
in order to compute the electric quadrupole moments of the nuclides !?!Sb and !?35b. 


Experimental Part 


The Sb I lines most conveniently analyzed lie in the near ultraviolet. These correspond 
to transitions from levels of configuration 5p?6s*Pyję, 32 to those of the ground configuration 
9p? ?P? 2, 3,2. The hfs of these lines was investigated by Badami (1932), who, using an antimony 
arc as light source, was able to assess the magnetic splitting constants for these levels. Basing 
on the data in Badami's paper, the present author decided upon the line A 3030A (5p? 2 12 > 
—>5p? 6s ?P„ 9) for measurements of the quadrupole coupling constant. 


FB h ; 
to spectrograp gam -— guartz window 


| to vacuum 


| - Iron 
| cathode 
| Ø 29/28 


liquid 


H u AL 
Ve) 
`+ 

Sb specimen 


s 


Fig. 1. Schüler type hollow cathode discharge tube used for exciting Sb I line (dimensions in mm). 


Arc excitation of the Sb spectrum was carried out by means of a Schüler type hollow 
cathode discharge tube connected to a vacuum apparatus. The inert carrier gas in the tube 
was purified in a circulating system. Å cross section through the tube and details of its con- 
struction are shown in Fig. 1. The specimen of metallic Sb in the form of a hollow cylinder 
of 45 mm length and an inner diameter of 5 mm was placed within the iron cathode coloed 
with water or liquid nitrogen. Using helium with a small admixture of argon for filling the 
tube, Sb T lines of high intensity are obtained. 


FTF 


The high resolving power instrument consisted of a Fabry-Perot interferometer in the 
form of a plane-parallel plate of fused silica with PbCl, and MgF, multilayer dielectric 
layers. The advantage of using PbCl, which is a compound possessing a high refractive 
index, for preparing the dielectric mirrors destined for use in the near ultraviolet, as well 
as the optical properties of such mirrors, have been investigated by Penselin and Steudel 
(1955). The present author prepared dielectric mirrors of PbCl, — MgF, in a vacuum appa- 
ratus wherein the thickness of the dielectric layer was measured during evaporation by a simple 
optical method?. 

By using nine-layer mirrors with a maximum of the reflection coefficient comprising 
the Sb line investigated, the resolving power of the interferometer necessary for exact hfs 
analysis was obtained; at the same time, the luminosity of the interferometer was several 
times better than that of aluminium mirrors. 

The interferometer was mounted in a parallel light beam in front of an 152-28 quartz 
spectrograph of medium dispersion. By using an interferometer in the form of a single 
plane-parallel plate, much tedious adjustment could be avoided, as otherwise dielectric 
mirrors having a maximum outside the visible range would render adjustment of the instru- 
ment excessively difficult. The interference fringes were mapped on the slit of the spectro- 
graph by means of a quartz-fluorite achromatic lens of f=25 cm. 

The distances between the hfs components of the A 3030A line were measured by photo- 
metrizing the contour of the structure with an MF-2 microphotometer. In order to obtain 
better accuracy in measuring the distances between the levels of hfs splitting of the term 
Ps, interference plates of various thickness were employed (3.65 mm and 6.00 mm). 
Moreover, the line 4 3267À (5p? ? 1/2 >5p? 65 ?P,5) having the level ?Pi in common 
with the line investigated was measured. The results for the 121Sb isotope are shown in Fig.2. 
The left-hand side contains the distances (in mK) as resulting from the interval rule, whereas 
the right-hand side gives the experimental mean values from many measurements. In Fig. 3, 
the energy levels of hyperfine splitting of the 2 3030A line are shown schematically, together 
with its hfs, without taking account of quadrupole interaction; the letters a, b...f denote 


8 
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5p? 65 Py 3 ——-t--- 
144.5 437 £ 1.5 
ہے ت ا‎ __ı__- 
96.4 
d 222 
theor. exper 


Fig. 2. Distances between the energy levels of hyperfine splitting of the term MM of 121 Sp, as computed from the 
interval rule, together with the experimental values (in mK). 


2 The simple method worked out by the present author for checking the thickness of dielectric mirrors to 


be used in near ultraviolet investigation will be published shortly. 
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Vig. 3. Energy level diagram and hyperfine structure of 3030 A line. Distances in mK. 


components of the isotope 121Sb, and A, B...F — those of 123Sp, It is to be regretted that 


there was no possibility of carrying out measurements of the distances between the hyperfine 
splitting levels in 123Sb with sufficient accuracy when using specimens of natural isotope 
composition, 
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For computing the quadrupole coupling constant of the term P32 of 1?1Sb, the energy 


(1) 


of the various levels can be written in the form (Kopfermann 1956) 


We Wy AĆ pp ACH - M+ TT 


2 2I(2I — 1)J(2J — 1) 
wherein 4 0 is the hfs splitting constant, 


B = eQ ووب‎ (0) — the quadrupole coupling constant, 


and C=FF+ 1) — I(I +1) — JJ 4-1. 


The effect of quadrupole interaction on the position of the energy levels for the term 


*Pyjg of the isotope 1?1Sb is shown in Fig. 4, wherein the continuous lines denote the positions 


> 


5p°6sP, 3 
2 
121 Sp (T= á) 


Fig. 4. Hyperfine splitting of energy levels for the term ركلة‎ and /—5/2 (12! Sb), account being taken of quadrupole 


interaction (continuous line for B>0, dashed line for B<0). 


of the levels for B>O, and the dashed lines are those of B<O. On comparing Figs. 2 and 4 
it will be seen that, in the case under consideration, the constant B is positive and the values 


of the constants 4,4, and By can be computed from the following equations: 


34 — £ B = (137 + 1.5) mK 
44 + 3 B = (206 + 1.5) mK 
Hence, for the term 5p? 6s ?P4, of the isotope 12!Sb, we have: 
Ajo, = (48.2 + 1.0) mK 
Bis = (16.7 + 2.5) mK 
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Computation of Q (P1Sb) and Q (Sb) 


With the constant B defined above (Eq. (1)), the quadrupole moment of the nucleus 
can be obtained from the following equation (Casimir 1936): 


QM E | 2 
eż); [r-3 - (3cos?0 — D] 
wherein the sum in the denominator extends over all electrons contributing to the constant B: 
in the present case, the sum extends over the three electrons forming the configuration p?s. 
Assuming strict Russel-Saunders coupling, and expressing rè by the spin-orbit interac- 
tion constant £j (Condon, Shortley 1957, Trees 1953), Eq. (2) can be rewritten as follows: 
BZ*H 
Ga 218cos*0 — 1) 
Here, B is expressed in 10-3 cm}, £,— in cm-?, Z* denotes the effective charge of the nucleus 
with respect to the outer electrons, and H=H (l, Z) is the relativistic correction to ¢). 


Q = x 0.7455 x 107% cm?. (3) 


In order to compute >] (3 cos? 0 — B we have to express the matrix element of the 
level p?s *Pzją as a whole by the matrix elements of the various electrons. To this purpose, 
applying the vector-coupling formulas (Condon, Shortley 1957, p. 228), we represent the 
eigenfunction of the level ?P,, as expressed in the representation (SŁMSM,), by means 
of the one-electron antisymmetric eigenfunctions in the representation (sl m, m), as follows?: 


V (pis PERAPY 1) = 76 [20(1*0*07) — (17005) ع‎ P(170*07)] (4) 


In all cases, the first two symbols in the brackets are for the two p electrons, whilst the third 
is that of the s electron. 
By Eq. (4), the matrix element of the level ?P,, can be written as follows: 
2 "a R 
( Ps), a 3 cos? — 1)|?Ps), ارم‎ ë [9 (L+|3 cos? — 1|14) + 
+ 5(0*|3 cos 29 — 110 + (0713 cos? — 107) + (1718 cos?ð — il. (5) 
Herein, the symbols of all functions refer to the p electrons, as the s electrons do not contribute 
to the matrix element (because of the spherically symmetrical charge distribution). 
The matrix elements (3 cos? 2— 1) for the p electron are now computed by making use 


of the angular part of the eigenfunction, ©,(l, mj), account being taken of the spin of the 
electron. We have now 


(Ps, m cos? — | pij, *) = (py, —"13 cost) — ll ps), =) = — i R 
q^) — - 1⁄2 a M __ 2 4 
(Pay, "213 cos?) — | Ps) = (Px a Ü lp; °) = 5 A 
(Dy, %|3 cos?) — 1 py”) => (py, Ká - Lips: 5%) = 0 
nią 2 1⁄2 VIAN ur -y 2/28 
(Pr, *|3 cos*à — 1| يرم‎ ^) = — (psy, ^|3 cos?0 — 1|py ^) = ب‎ 2s MEE (6) 


£ aj 81% . . . H 0 . . 
R, and S, are relativistic corrections according to the notation of Casimir (1936). 


? The notation is that of Condon and Shortley (1957). 
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In order to be able to apply the foregoing elements in Eq. (5), they should be expressed 
in the representation (sl m, mj): 
(I*3cos?0 — 1|1+) = — > R; 
+ 2 — áp 
(07|3 cos? — 1|0*) = R 
(0-13 cosż8 — 107) = 4 (R; + 25,) 
| (1-|3 cos? — 1|1-) = $ (R, — 45). (7) 
On substituting (7) in Eq. (5), we have ultimately , 
EP 1 (3 costó — DI. = 2 (— به‎ +108,) 
For antimony (Z=51), R, = 1.1182, S,= 1.2015 (Kopfermann 1956), whence the numerical 
value of the matrix element under consideration is 
Á CA] D (3 cos? BRA DI = 0.4843 
On substituting in Eq. (3) the constant B=(16.7 +2.5) mK as found for the 5p? 6s "Psia 
level of the isotope 121Sb with Z* = Z—4=47 for the p electrons, H=1.0657, 6, =2400.0 cm-1 
and 2j (3 cos28—1) = 0.4843, we have 
Q(-15b) = (— 0.54 + 0.08) 10-24 cm? 
The electric quadrupole moment of the isotope !?3Sb can be computed from the forego- 
ing value of Q,,, and the ratio Qg1/Q123, which is known to a high degree of accuracy from 
resonance measurements. According to Wang (1955), Qis4/Q45, — 1.274745 +0.000010, 


whence 


Q(1238b) = (— 0.69 + 0.10) 10-24 cm? 


Sternheimer's (1950, 1951, 1952) corrections have been omitted. 


Quadrupole moment values for 1*!Sb and !??Sb as hitherto obtained by various authors 


The values hitherto obtained by various authors for the quadrupole moments of either 


stabile isotope of antimony are assembled in the following Table: 


121 Sh 123 Sh ae Á des 
0x1024 cm? Qx102 cm? Method of measuremen uthor 
0 40 from hfs of Sb II lines Tomboulian, Bacher (1940) 
—0.3--0.2 —1.2-0.2 from hfs of Sb II lines Murakawa, Suwa (1949) 
—0.35 to —0.5 —0.45 to —0.62 from microwave measurements | Loomis, Strandberg (1951) 
— 1.3 — 1.7 from hfs of Sb II lines Sprague, Tomboulian (1953a, 1953b) 
—0.53--0.10 —0.68+0.10 from hfs of Sb II lines Murakawa (1954, 1955) 
— 0.71 to —14 —0.91 to —1.8 from microwave measurements | Jache eż al. (1955) 
(—0.8) (—1.0) 
—0.54--0.08 —0.69+-0.10 from hfs of Sb I lines present author 
—0.68 — 1.0 theoretical values Wageningen, Boer (1952) 


— 0.26 — 0.39 theoretical values Horie, Arima (1955) 
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From the foregoing Table, the present author's results are seen to be in almost exact 
coincidence with those of Murakawa, notwithstanding the fact of their having been obtained 
from the arc (and not spark) spectrum of Sb. This should seem to prove that the results of 
Murakawa and of the present author are correct. 

Finally, the author wishes to thank Professor H. Niewodniczański, Director of the 
Institute of Physics of the Jagellonian University, for his valuable advice and discussions 
throughout the present investigation. The author also thanks Docent T. Skaliński, Director 
of the Institute of Experimental Physics in Warsaw, for kindly providing him with the 
interference quartz plate of thickness -3.65 mm, which greatly facilitated the measurements. 
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LETTERS [O THE EDITOR 


SPIN-SPIN RELAXATION TIME IN A MODEL OF NON-UNIFORM 
FERROMAGNETIC 


By JANUSZ MORKOWSKI 
Department of Ferromagnetics, Institute of Physics of the Polish Academy of Sciences, Poznań 


(Received April 12, 1961) 


A brief report on spin-spin relaxation time calculated for a simplified model of non- 
-uniform ferromagnetic is presented here. Details of the caleulation and some other results 
will be published in future. 

The model of non-uniform ferromagnetic considered here is an array of twod ifferent 
spins, S and S’, randomly distributed over lattice sites with relative concentrations c and 
1 —c, respectively. The lattice is assumed crystallographically perfect, with a cubic symmetry. ' 

Magnetic properties of the model are determined by its Zeeman energy in an applied 
constant magnetic field and by mutual interactions between spins, i.e. exchange interactions 
and anisotropic ones. Exchange interactions between two neighbouring spins have the 
usual form of the scalar product of their spin operators with coupling constant (exchange 
integral) equal J for two nearest spins S, J’ for the interaction of spin S with neighbouring 
spin S’ and finally J” for interaction of two nearest spins S’. The anisotropic magnetic dipolar 
interactions as well as the short range pseudodipolar ones are taken into account. The pseudo- 
dipolar interactions are confined to the nearest neighbouring spins only, with coupling 
constants D, D” and D' for interactions between two spins S, two spins S’ and for S with 54 
respectively. 

The spin-spin relaxation time is calculated in the spin-wave approximation. Transverse 
components of the spin operators are replaced by harmonic oscillator operators according 
to the rules 5; — )25( ور‎ Sit > (28')% nf — a generalization of the ones employed 
for an uniform system (cf. Morkowski 1960, Eq. (3.18)). As usually, operators «y, 1; 
are expanded into their Fourier components. The diagonal part of the Hamiltonian 
of the spin system, averaged over all distributions of spins S, S' determines the 
energy spectrum of the spin waves. In a finite sample (spherical shape of the crystal is 
assumed here) magnetic dipolar interactions lead to the known degeneracy of the spin-wave 
(Clogston et al. 1956). Interaction part of the Hamiltonian causes transitions 


(783) 


spectrum 
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between various spin-wave states determining relaxation time of the system. We discuss 
here only the scattering processes of the lowest order, i.e. two-spin-wave processes governed 


by the term >) Wy; y, of the interaction Hamiltonian (y7, y; are the spin wave creation and 


4,0 : > 5 
annihilation operators). These processes, which do not conserve momenta of interacting spin 


waves, are forbidden in all systems invariant for translations, for instance in the uniform 
system. In disordered ferrites these two magnon processes lead to an inverse relaxation 
time comparable with observed resonance linewidth (Clogston et al. 1956, Callen and Pittelli 
1960). 

The spin-spin relaxation time is calculated in the usual way, employing the kinetic 
equation for small deviation of the spin-waves system from the thermodynamic equilibrium. 
The relaxation time for long spin waves u (with wave vector magnitude u equal to the in- 
verse skin-depth), excited by microwave magnetic field in the ferromagnetic resonance is 
calculated. The results depend on a direction of m. Therefore, for definiteness, the direction 
of perpendicular to the applied magnetic field (directed along easy axis of magnetization) 
was chosen. Simple final formulae (valid for large magnetic field H) are given below for two 
special cases: 

a) In the case of small magnetic impurity (1—c<l), with the assumption J—J' =J”, 
D=D'=D', the inverse relaxation time is 


1“ 12 S'% |* (drupM,JS)% 5 
T = 24n c(l c) (1 za h 1 STAŃ $ (1) 
for 
daupMy 
x ; zul 2 
%= 2ugH + c(1—c) (5⁄2 — S 2)? Jz a (4) 


where z is the number of nearest neighbours, M, — the saturation magnetization, n is equal 
1, 2 and 4 for simple, body centered and face centered cubic lattices, respectively. 

b) Relaxation time for the case of dilute ferromagnetic (for high dilution, cl, the 
formula below is not applicable) 5'— 0, J=J"=D=D'=0, is 


a 1.102 S\% 
1-28 1-9: (^) (Arun M, JES) (1 e E 


Ann h 
1 Agu My \* (4eug M, Jees) 9 
2, K 3 0 0 € 
+ Te (1— c) ( TES ) ^ 1 8 6+... (3) 
for 
2n Moc 
K= Hn (4) 


d is equal 2, 0, 1 respectively for simple, body centered and face centered cubic lattices. 
The quantity Jc? plays a role of an “effective” exchange integral, which can be estimated 
from experimental value of Bloch 732 law constant Cz= (0.0587/nS) (k/Jc2S)3/2, 

The proper application of the results above is to disordered ferromagnetic alloys. We 
can, however, tentatively apply the model considered to ordinary ferromagnetic metals, 
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such as for instance nickel. As a crude approximation to the real situation we consider the 
model of a nickel crystal, in which Ne (c=0.606=effective magneton number per atom) 
atoms have spin S=1/2 and remaining N(1—c) atoms — spin S'=0. For values for J, D 
estimated from Bloch 75/2 law and from anisotropy constant the value of 1/7 is of the order 
of magnitude of measured resonance linewidth. It is, however, difficult to draw out more 
definite conclusions as the spin-spin relaxation processes are not the sole contribution to 
the experimental resonance linewidth. 


Sincere thanks are due to Professor S. Szezeniowski for his interest and continuous 
encouragement during this work. 
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ON THE QUANTUM THEORY OF SPIN-WAVE RESONANCE IN A FINITE 
FERROMAGNETIC LINEAR CHAIN 


By A. R. FERCHMIN 
A. Mickiewicz University, Institute of Theoretical Physics, Poznań 
(Received April 27, 1961) 


In this letter there will be presented a method of diagonalization of spin Hamiltonian of a ferromagnetic 
inear chain taking into consideration instead of the usual cyclic boundary condition a more realistic one. As an 


application standing spin-wave resonance (Kittel 1958) will be considered. 


I. Diagonalization of the Hamiltonian 


We restrict ourselves to a one-dimensional line of N spins placed on the z-axis. A mag- 
netic field H(h,, h,, Ho) is applied. H, denotes the static magnetic field in the z direction, 
h, and h, are components of the uniform rf magnetic field h of amplitude ho, circularly | 
polarized. The Hamiltonian of the system consists of Zeeman energy of spins in the applied 
field, of the mutual interactions between spins 1. e. exchange interactions with coupling 
constant (exchange integral) J, and anisotropic short range pseudodipolar interactions with 
coupling constant D,, (see, e. g., Morkowski 1960). Only nearest neighbour interactions are 
considered. The exchange integral near the “surface” Jg =Jn-ın =J FAJ, Jiri =J a; SJ 
for all others site numbers / (cf. Bontch-Bruevitch 1956). Similarly, we assume that D} = 
=Dy 1 N"D+4D, D,;,,=D;_,,=D for all others 2. According to Dyson’s (1956a, b) 
calculations we introduce “physical” and “ideal” states using localized oscillator amplitudes 
ررد‎ and n. The non-Hermitian Hamiltonian can be expressed in terms of those operators. 
For details see Morkowski (1960). The bilinear part of the Hamiltonian can be diagonalized 
by means of the following Fourier transformations: 

NEC > a =E 5) ni. (c — normalizing factor) (1) 
k=1 k=1 
The (real) functions f, satisfy the following set of difference equations: 
fall 0) 4 flw E oO ôm)] tfpa (1 ôm) = 0; 1, «as, N (2) 


where: 
_ 9 J Deon = e/5 pnp per AF AD 
J-D $ J-D 
(787) 
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lg — Bohr magneton, S — magnitude of spin vector, € — energy level to be determined 
from the secular equations (5) by means of relation (4). 
One finds from (2) (cf. Ferchmin 1960): 


fı = A [cos 10 + tan 4 (N — 1) © sin 10] (for symmetric modes), 
fui = A [cos 10 — cot } (N — 1) © sin 1©[ (for antisymmetric modes), 
i+1=1, .., N. (3) 
The parameter © is defined by the relation 
— 2 cos © =w (4) 
The secular equations have the following form: 
cos 4 (N+ 1( © — (1 — o) cos + (N — 1) 0 = 0 (for symmetric modes), 
(sin ©)-1 [sin 1 (V+ 1) © — (1— o) sin HN —1)0]=0 (for antisymmetric modes) (5) 


` The parameter © is identical with the wave number of a standing spin wave multiplied 
by the lattice constant. 


II. Spin-wave resonance 


The part of the diagonalized Hamiltonian containing the external rf field is: 


N N 
Hy = OC 2; nt DSi (6) 
=1 ia 


The derivation of transition probability for non-Hermitian Hamiltonian and for ideal states 
with indefinite metric by means of perturbation theory was given by Morkowski (1960). 
Using his methods we get the resonance transition probability per unit of time from the 
ideal state |d) to the state |a) caused by rf field: 


N N | 
1 2 
Via = "t EA 25 (az + 15 2; tr ) Öle — eg - ho) (7) 


t, t, energy levels of states |a) and |d), respectively, 
a, — number of spin waves with wave number k in |a)-state. 
The sum MY, ff in (7) equals zero for antisymmetric spin-wave modes, and 
~ in 3 NO, ` 1 = 

" sin + k oy Jf v 

D>A— 2 = A — sin 0 8 
Daf sin 4 Ox cos 1 (N — 1) O, 2 | 2 | (8 
for symmetric spin-wave modes. 

We conclude, that spin waves with 9,70, e. g., for k#0, can be excited when o (depending 
on both isotropic part of the surface energy and the anisotropic one) #0. Moreover, the 
so-called oscillator strength is proportional to k“? for small wave numbers k. The latter result 
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disagrees with that obtained by Kittel (1958) semiclassically and by Pincus (1960) using 
Bloch spin wave theory. These Authors have found the oscillator strength proportional 
to k-! (provided that k is small). 

| Calculations referred above can be extended to cover the two- and three-dimensional 
cases. A detailed investigation of this point is currently under way. 

The author is very indebted to Professor 5. Szczeniowski for his interest and for reading 
the manuscript. Thanks will be due to Dr. J. Morkowski and to Dr. L. Kowalewski for 
helpful discussions. 
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REVIEWS OF BOOKS 
“THE ENCYCLOPEDIA OF SPECTROSCOPY" 


GEORGE L. CLARK 
Reinhold Publishing Corporation New York Chapman and Hall, Ltd., London, 1960, pp. 787, ... 


“The Encyclopedia of Spectroscopy” edited by G. L. Clark contains contributions by 
121 eminent specialists in the various branches of spectroscopy. 

“The Encyclopedia of Spectroscopy” comprises the following Sections in alphabetical 
order: Absorption Spectroscopy (Spectrophotometry) — visible and ultraviolet, Band 
Spectroscopy, Beta-Ray Spectroscopy, Differential Thermal Analysis, Electron Paramagnetic 
Resonance Spectrometry, Emission Spectroscopy — Light, Flame Photometry, Fluoro- 
photometry and Phosphorimetry, Gamma-Ray Spectrometry, Infrared Emission Spectro- 
scopy, Infrared Spectrophotometry, Mass Spectrometry, Microwave Spectroscopy, Mono- 
chromators, Neutron Spectrometry, Nuclear Magnetic Resonance Spectra, Raman Spectro- 
scopy, Solar Spectroscopy, Vacuum Spectroscopy, X-Ray and Gamma-Ray Absorption , 
Photometry (Absorptiometry), X-Ray Characteristic Absorption Spectrometry, X-Ray 
Diffraction or Crystal Spectra, X-Ray Emission Spectrometry. Within the various Sections, 
the items are arranged in alphabetical order. 

The items are conceived so as to give a brief review of the historical development of 
the branch of spectroscopy under consideration, some quite general theoretical treatment, 
a rather comprehensive critical account of the most recent experimental methods with 
data on their degree of accuracy and limits of applicability, as well as a review of the appa- 
ratus in use, together with pictures of the most commonly used factory-made devices. In 
addition, each Section brings a detailed discussion of the scientific, technical or medical 
applications of the respective branch of spectroscopy. The highly comprehensive list of 
references attached to each problem is especially valuable, as it contains the most recent 
papers published, and even a number of unpublished items. 

“The Encyclopedia of Spectroscopy” lays especial stress on practical problems related 
to the applications of absorption spectroscopy, emission spectroscopy, infrared spectroscopy 
and X-ray spectroscopy. Thus, “The Encyclopedia of Spectroscopy” can be of great use 
as a work of reference in metallurgical, mineralogical, biological, medical, geological, agro- 
nomical, astrophysical and other laboratories. The vast amount of references it contains 
makes it a highly desirable help in physical or chemical spectroscopic research work. 


Maria Danuta Kunisz 
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ERRATA 


„Ihe Critical Problems for Multilayer Slab Systems“, A. Kuszell (Acta Physica Polonica, 20, 567 (1961)). 


it should be 


tanh 1 
tanh“! 
V, (x, u) = V,(— x. — p) 
N,(v.) 
x0): ZX) 
B3(— v) x (— ») 
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